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ABSTRACT 


The  intensities  of  the  second  harmonic  light  in  total  reflection 
and  in  transmission  from  uniaxial  KDP  cyrstals  immersed  in  a«n  optically 
denser  liquid  1- bromonaphthalene  have  been  measured  as  a  function  of 
angle  of  incidence  of  the  fundamental  beam  from  a  mode  locked  Nd:glass 
laser.  In  the  experiment ,  va rious  crys tallographic  orientations  of  the 
KDP  crystal  were  used.  The  results  agreed  well  with  the  theoretical 
predictions  of  Bloembergen  and  Pershan  .  In  particular,  the  existence 
of  a  nonlinear  Brewster'  s  angle  of  a  transparent  medium  (KDP),  with 
nonlinear  polarization  in  the  plane  of  reflection,  has  been  first  demon¬ 
strated.  The  transmitted  second  harmonic  intensity  under  phase  matched 
condition  by  birefringence^  ’  \  in  the  oblique  direction  was  found  to 

vary  more  than  ten  orders  of  magnitude  in  the  angular  range  of  interest. 

The  transmitted  second  harmonic  intensity  under  noncollinear 

phase  matched  condition  by  means  of  two  beam  spatial  mixing  (TBSM) 

in  the  KDP  crystal,  using  Q-switched  Ndrglass  laser  pulses,  was 

observed  as  a  function  of  the  angle  of  incidence.  The  result  was  in  good 

(12  ) 

agreement  with  the  theory  of  Bloembergen  and  Pershan  .  As  a 
consequence,  this  result  was  utilized  to  measure  the  second  order 
intensity  auto-correlation  function  of  Nd:glass  laser  radiation.  The 
picosecond  pulsewidth  of  Nd:glass  laser  radiation  was  measured  to  be 
5.29  picosecond. 
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CHAPTER  I 
INTRODUCTION 

A.  General  Introduction 

Since  the  advent  of  LASER  (1)  (Light  Amplification  by  Stimulated 

Emission  of  Radiation)  or  optical  maser,  the  coherent,  high  intensity 

and  monochromatic  light  sour ces  are  available  for  nonlinear  optical 

processes  in  materials.  Franken  and  his  co-worker ,2  >.  in  1,61,  performed 

an  artfully  devised  experiment  in  which  they  proved  that  the  ruby  laser 

light  of  wavelength  6943A  focused  within  a  quartz  crystal  contained  on 

emerging  a  small  admixture  of  a  second  harmonic  of  a  wavelength 
o 

3471.  5A.  The  discovery  of  this  phenomenon,  which  has  its  origin  in  the 

nonlinear  optical  properties  of  matter,  became  the  starting  point  of  vast 

development  further  enhanced  by  the  rapid  evolution  of  laser  technique. 

This  has  led  to  the  detection  of  a  number  of  novel  optical  effects, 

which  have  since  been  discussed  in  various  monographs^3,  ^  and 
,  (5-10) 

review  articles  .  Furthermore,  the  theoretical  essentials  of 

coherent  nonlinear  processes  are  to  be  found  in  the  fundamental  work  of 
Armstrong  etal(U),  Bloembsrgen  and  Pershan'12',  Akhamonov  and 

/I  o  \ 

Khokhlov  and  certain  other  authors ^ 

When  low  intensity  light  propagates  in  a  transparent  medium,  it 
does  not  affect  the  optical  properties  of  the  medium  because  the 

electric  polarization  P(r,t)  induced  in  it  is  a  linear  function  of  the 

i  •  .  *  -* 

electric  field  E(r,  t)  of  the  light  wave.  A  linear  relationship  between 

-*  -*  -♦  ■+ 

P(r,t)  and  E(r,t)  results  immediately  in  I^rentz'  s  classical  electron 
(18) 

theory  .  To  a  linear  approximation,  the  solutions  of  Maxwell'  s 
equations  satisfy  the  principle  of  superposition,  which  states  that  the 


1 


electromagnetic 


waves  simultaneously  transversing  a  linear  medium 


? 


propagate  independently  without  interacting  upon  one  another.  In  linear, 
optically  transparent  medium,  electromagnetic  waves  accordingly 
propagate  without  distortion,  the  refractive  index  n  depending  solely 
on  the  properties  and  thermodynamic  state  of  the  medium  but  not 
depending  on  the  incident  light  intensity.  All  optical  phenomena  depend 
linearly  on  the  electric  field  strength  as  long  as  the  medium  under 
investigation  is  probed  by  low  intensity  light  wave  and  is  not  acted  upon 
simultaneously  by  external  field. 

However,  when  a  transparent  medium  is  transversed  by  light  of 

very  high  intensity,  i.  e.  laser,  its  refractive  index  n  and  electric 

permitivity  become  dependent  on  the  intensity  of  the  incident  light. 

This  has  indeed  been  directly  observed,  with  the  help  of  giant  pulsed 

ruby  lasers,  by  various  authors  (2  ’  19_  21  Lasers  operating  continuously 

or  pulsed  radiate  a  parallel  and  monochromatic  light,  coherent  in  space 

and  time,  and  contain  a  flux  of  immense  energy  density.  The  very  high 

light  intensities  of  the  beams,  never  achieved  by  ordinary  light,  cause 

the  index  of  refraction  to  depend  on  the  electric  field  E(r,t)  of  the 

beam.  An  electromagnetic  wave  of  such  strength,  inducing  optical 

nonlinearity  in  the  medium,  itself  undergoes  distortion  when  propagating 

through  the  medium.  By  Fourier  spectral  analysis,  the  original  wave 

of  fundamental  frequency  uj  is  now  additionally  accompanied  by  harmonic 

components  of  double  frequency  2  14  triple  frequency  3  U.'  and  higher 

(3-6,  11-13) 


harmonics . 


.  Furthermore,  due  to  the  availability  of  high 


intensity  and  ultrashort  light  pulse  s,  various  nonlinear  optical  processes, 
e.g.  self  focusing  phenomena,  self  induced  transparency,  stimulated 
Raman  emission,  photon  echo,  have  been  observed  and  well  described 
in  various  articles^2  24)  and  reviewed  article09,  10,  17).  the  present 

work,  which  will  be  presented  here  thereafter,  only  second  harmonic 
generation  in  reflection  and  transmission  from  an  uniaxial  crystal  KDP 
using  Nd:glas s  laser  radiation  as  fundamental  beam,  will  be  described  in 
detail,  and  to  a  lesser  extense,  an  attempt  to  generalized  the  effect  to 
third  order  is  also  included. 

B‘  Review  of  Previous  Experiments 

—  rm°nic  Generation  in  Reflection.  Theory  of  interactions 
between  light  waves  in  a  nonlinear  media  was  given  by  Armstrong. 
Bloembergen,  Ducuing  and  Pershan'11’.  Consequently,  the  theory  of 
light  wave  at  the  boundary  of  nonlinear  media  was  given  by  Bloembergen 
and  Pershan  .  Second  harmonic  generation  in  reflection  was  first 
observed  by  Ducuing  and  Bloembergen  in  1963<25’  26).  Both  the  directional 
and  the  polarisation  properties  of  second  harmonic  reflection  from  III-V 
semiconductor  mirrors,  e.g,  GaAs,  have  been  verified.  Furthermore, 
Bloembergen  has  predicted'12  ’  and  demonstrated'26’  SHC  in  reflection 
by  using  two  fundamental  beams  of  Q-switched  ruby  laser  spatial  mixing 
in  GaAs  mirror.  Chang  and  Bloembergen'27’  have  verified  the  laws  for  the 
reflected  intensity  of  second  harmonic  light,  hi  the  experimental 
verification  they  used  various  crystallographic  cuts  of  cubic  crystal 
GaAs  mirrors  upon  which  the  fundamental  beam  of  Q-switched  ruby  laser 


was  incident, In  that  experiment  Chang  and  Bloembergen  have  first 

demonstrated  the  existence  of  Nonlinear  Brewster  Angle  for  the 

absorbing  nonlinear  medium.  Bloembergen  and  Lee,28>  observed  the 

internal  reflected  second  harmonic  intensity  (SHI)  generated  at  4860A 

by  an  incident  beam  totally  reflected  by  nonlinear  medium  NaCIO  and 

3 

KH2  P°4  (KDP>-  immersed  in  an  optically  denser  linear  fluid  (1- 
bromonaphthalene ).  The  incident  beam  was  the  stimulated  Stoke  beam 
induced  by  Q-switched  ruby  laser  in  H2  gas.  In  accordance  with  theory 

given  by  Bloembergen  and  Pershan^12^ 

S  ana  ^ersban  ,  this  experiment  shows  anomalously 
high  reflected  harmonic  intensity  because  of  the  enhancement  of  the 

intensity  by  phase  matching  at  the  critical  angle  of  incidence  for  total 
reflection.  Furthermore,  a  detailed  discussion  and  experimental 
verification  to  the  theory'12 »  concerning  total  reflection  and  two  beam 
spatial  mixing  on  NaClQ,  has  been  given  by  Bloembergen,  Simon  and 

t _ (29) 


>ergen  and  Lec^  they 


The  first  experiment  in  obtaining  an  ultraviolet  third  harmonic 
generation  in  reflection  was  done  by  Bey,  Giuliani  and  H.  Rabin(30). 
With  the  help  of  previous  work  of  Bloembergen  and  Lee(28),  they 
introduced  in  the  experiment  anomalous  dispersion  into  normally 
unmatched  medium.  They  obtained  phase  matched  third  harmonic 
generation  (THG)  at  3530A  of  the  Nd:glass  laser  line  by  introducing  dye 
molecule  of  fuchsin  red  into  the  liquid  medium  of  hexafluoroacetone 
sesquihydrate.  Wang  and  Baarden !  1  reported  studies  of  THG  reflected 

from  the  boundary  of  solids  including  metal,  semiconductors,  alkali 


halides  and  glass.  Recently,  Burns  and  Bloembergen(32 )  used  pico¬ 
second  pulse  from  Nd:glass  laser  to  study  reflected  THG  in  absorbing 
media  of  cubic  or  isotropic  symmetry. 

^ ■  Harmonic  Generation  in  Transmission.  The  first 

experiment  of  second  harmonic  generation  was  performed  by  Franken 

(2) 

et  al  in  the  transmitted  direction.  They  used  ruby  laser  beam  focused 
in  a  piece  of  quartz  and  observed,  after  separation  the  fundamental  by  a 
prism,  a  second  harmonic  of  wavelength  3471.  5A.  Since  then,  the 
second  harmonic  of  light  has  been  observed  in  various  piezoelectric 
and  other  nonlinear  material.  Giordmaine(33),  Maker  et  al(34)  in  1962, 
demonstrated  independently  the  use  of  birefringence  of  uniaxial  crystal 
to  achieve  phase  matching  condition  under  which  the  second  harmonic  of 
light  is  enhanced.  Armstrong  et  al(11)  and  Kleinmann(14)  introduced 
the  concept  of  coherent  length  along  which  the  fundamental  and  second 
harmonic  beams  will  travel  in  step  with  respect  to  each  other.  Second 
harmonic  generation  with  double  refraction  has  been  given  in  detailed 
by  Boyd  et  al(35).  Miller(36)  used  Nd:CaWG4  laser  to  observe  second 
harmonic  of  light  from  BaTi^,  KDP,  ADP  CdS  and  strongly  nonlinear 
LiNbG3  and  LiTa  .  Geusic  et  al(38)  observed  SHG  at  5320A  in 

Ba2NaNbOi5  and  K6Li4Nb°3  using  slightly  focused,  continuously  pumped, 
repetitively  Q-switched  NdrYAG  laser.  Using  ruby  laser  and  Ndrglass 
laser  Kurtz  et  al^  ^  and  Nath^  ^  made  a  novel  study  of  SHG  from 
strongly  nonlinear  materials  HI03  and  LilO  respectively.  Askin  et  al(41) 
observed  continuous  SHG  using  the  infrared  transition  1. 1526p  of  the  He-Ne 


gas  laser  with  a  focused  and 


a  non-focused  beam  to  measure  the  second 


order  susceptibility  component  d ^  of  KDP.  With  the  help  of  temperature 
tuning,  phase  matched  SHG  in  LiNbC^without  double  refraction  using 
infrared  transition  1.  1526p  of  He-Ne  gas  laser  was  demonstrated  by 
Miller  .  Bio embergen  et  all  demonstrated  two  beam  spatial 
mixing  in  NaClC^  near  critical  angle  of  incidence  and  observed  SHG  in 
transmitted  direction.  The  equation  for  the  nonlinear  phase  matched 
condition  which  requires  all  three  propagation  vectors  of  the  fundamental 
and  of  the  second  harmonic  have  to  form  a  closed  loop  was  first 
pointed  out  by  Tien(4i)  and  later  by  Boy.d(44).  Phase  matched  SHG  with¬ 
out  double  refraction  (or  noncollinear  phase  matched  SHG)  in  ADP, 

KDP  using  He-Ne  line  of  1.  152  6  p  was  demonstrated  by  A.  Askin  et  al(45). 
Later  noncollinear  phase  matched  SHG  in  Lil03  was  performed  with 

Q-switched  Nd:YAG  laser  by  Shinsuke  Umegaki  et  al(46).  Recently 
(47) 

Aggarwal  achieved  noncollinear  phase  matching  in  GaAs  by  using  a 
carbon  dioxide  laser. 

The  first  experiment  in  obtaining  an  ultraviolet  third  harmonic  at 
o 

2314  a  using  ruby  laser  is  due  to  Terhune  et  al(20)  in  calcite  and 
subsequently  in  cubic  crystal  LiF,  KC1,  CaF  and  liquids  by  Maker(48,  49). 
Bey  et  al(  )  provided  the  earliest  experimental  evidence  that  phase 
matching  can  be  achieved  in  harmonic  processes  by  the  introduction  of 
anomalous  dispersion  into  a  normally  unmatched  medium.  They  obtained 
phase-matched  THG  at  3530A  of  the  Nd:glass  laser  line  by  introducing 
dye  molecule  (fuchsin  red)  into  a  liquid  medium  (hexafluoroacetone 


s esquihydrate ).  Their  experiment  was  extended  by  Chang  and 
(57) 

Galbraith  to  other  solvents  of  different  index  mismatch  and  to 

another  dye,  methylene  blue,  Bey  et  al^  ^  reported  interesting  THG 

studies,  which  prove  that  a  linearly  polarized  laser  beam  generates  a 

third  harmonic  signal,  while  a  circularly  polarized  beam  does  not,  in 

accordance  with  theory.  Ward  and  New(53)  observed  THG  at  2314^  in 

gases  He,  Ne,  Ar  using  a  focused  ruby  laser  beam.  Harris  et  al(54) 

made  a  theoretical  analysis  and  proposed  THG  in  phase  matched  alkali 

metal  vapors,  e.  g.  Li,  K,  Na,  Rb,  Cs.  Cons  equently  Young (55  >, 

....  (56) 

Miles  demonstrated  THG  in  phase  matched  Rb  vapor  and  in  the  alkali 
metal  vapors  respectively. 

3'  Measurement  of  Picosecond  Pulsewidth  by  Means  of 
Nonlinear  Optical  Method.  The  Ndiglass  laser  has  recently  become 
an  important  research  tool.  Not  only  does  it  constitute  one  of  the  most 
powerful  sources  of  coherent  radiation  in  terms  of  peak  power,  but  it 
also  distinguishes  itself  as  a  source  of  ultrashort  light  pulse.  The 
emission  of  trains  of  short  pulses  produced  by  mode  locking  of  this 
laser  was  first  observed  by  DeMaria  et  al(56'  57).  Recent  review 
articles  of  mode  locked  laser  pulses  and  picosecond  laser  pulses  have 
been  given  by  DeMaria  (58),  DeMaria  et  al,59)  and  Smith(60).  fa  ruby 

i  (61  ) 

,  mode  locking  could  be  achieved  by  Q- switching  the  laser 
with  a  saturably  absorbing  dye.  The  duration  mode  locked  laser  pulses 
from  Ndiglass  laser  was  measured  by  Armstrong 162  1  by  means  of  SHC 
in  GaAs,  and  the  result  of  pulsewidth  measurement  was  about  4-6  pico- 


seconds  (ps).  The  introduction  of  another  technique  for  measurement  of 
the  pulsewidth,  so  called  two-photon  fluorescence  (TPF)  method,  was 
given  by  Giordmaine  et  al*  \  Rentzepis^64^.  In  the  TPF  method  care 
has  to  be  taken  in  connection  with  the  contrast  ratio  between  the  back¬ 
ground  and  the  overlapped  region  of  the  autocorrelation  picture. 

Various  theories  were  proposed  to  explain  the  picosecond  substructure 
observed  in  Nd:glass  laser.  The  simplest  proposal  made  independently 
by  several  authors (  69 \  was  that  the  picosecond  peaks  observed  in 

TPF  arose  from  the  short  duration  fluctuations  present  in  spontaneous 
emission.  The  theory  of  two-photon  fluorescence  and  another  nonlinear 
optical  coincidence  technique  (harmonic  generation  technique)  has  been 
treated  in  detail  by  several  authors,  particularly  rigorous  and  exhaustive 
treatment,  are  given  by  Rowe  and  Li(?0),  and  Picard  and  Schweitzer(71). 
The  difficulty  in  observating  the  contrast  ratio  3:1  was  removed,  and 
the  experimental  result  in  pulsewidth  measurement  in  agreement  to  the 
theory  was  performed,  by  means  of  TPF  with  very  thin  cell,  by  Shapiro 
and  Duguay  .  They  obtained  the  pulsewidth  of  Nd  rglass  laser  of  0.  4 
ps.  The  three  photon  fluoresence  technique  for  measurement  of  ultra- 
short  pulse  was  performed  by  P.  M.  Rentzepis  et  al(73),  and  the 
experimental  result,  in  agreement  to  the  theory,  of  the  contrast  ratio 
of  10:1  was  obtained.  The  exhaustive  review  of  study  of  the  Nd:glass 
Laser  Radiation  was  given  by  Duguay  et  al^74^. 

The  optical  harmonic  generation  technique,  which  gives  background 
free  in  autocorrelation  measurement  of  the  ultrashort  pulse,  was  first 
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performed  by  Armstrong'62 ».  *,  his  experiment  he  utiIized  ^  unjque 
polarization  property  of  of  GaAs,  given  by  Chang  et  al<27>,  for  SHG  in 

reflection.  Armstrong  obtained  the  pulsewidth  of  Nd:glass  laser 

radiation  of  4-6  ps.  Weber^^  „  ■  T,_.„  • 

’  using  KDP  in  experimental  set  up 

similar  to  famous  Michelson-Morley  experiments,  determined  the 
pulsewidth  of  Nd:glass  laser  radiation..  He  obtained  the  pulsewidth  of 
8  -12  ps.  Giordmaine  et  al  <76 '  introduced  a  new  technique  for  determining 
of  a  narrow  pulsewidth  by  using  two  beam  spatial  mixing  on  a  crystal  KDP 

th6ir  eXPerimen‘emPIOying  “  *"*«■•  light  burst  in  the  stimulated  Rama, 
effect  in  CS2  using  Q_ switched  ruby  laser,  the  backward  Raman- Stoke  wave  was 
equally  splitted  into  two  parts  with  a  proper  delay  time  T.  Then  the 
two  beams  were  recombined  in  KDP  crystal  oriented  to  allow  phase- 
unmatched  second  harmonic  generation.  The  average  Stoke  pulsewidth 
was  found  to  be  30  ps.  Treacy<77>,  using  a  grating  pair  for  pulse 

compression  technique,  obtained  the  pulsewidth  of  Nd, glass  laser  radiation 

Of  0.44  ps,  in  agreement  to  theoretical  limit  n.  • 

1  1116  unique  polarization 

Properties,  reported  by  Bey  et  «,<“>.  for  T„C  from  a  phase  matched 
solution  consisted  of  dye  fuchsin  dissolved  in  hexafluoroacetone 
sesquihydrate ,  was  promptly  utilized  by  Eckardt  and  Lee™  for  third 
order  autocorrelation  measurement  of  Ndtglass  laser.  The  result  of 
their  experiment  yields  the  pulsewidth  of  Ndtglass  laser  radiation  of  0.  7 
ps.  in  agreement  to  theoretical  limit.  Recently,  Jayaraman  and  Lee'79’  80> 
made  an  observation  of  two  photon  and  three  photon  conductivity  in  GaAs, 

CdS  and  CdS^  5-Se^  ?  using  Q.  switched  and  mode  locked  pulse  from 
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Ndrglass  laser,  and  demonstrated(8l)  the  potentiality  of  multiphoton 
conductivity  for  measuring  mode  locked  pulse  from  Nd:glass  laser. 

c-  Aim  of  the  Experiment 

It  is  the  purpose  in  this  work  to  study  and  verify  the  theoretical 
prediction  in  the  area  of  second  harmonic  generation  given  by 
Bloer  ibergen  and  Pershan'12’.  Furthermore,  the  experimental  results 
concerning  unique  polarization  properties  and  phase  matched  condition 
of  uniaxial  Potassium  di  Hydrogen  Phosphate  (KDP)  crystal  will  be 
utilized  for  measurement  of  picosecond  pulsewidth  of  Nd:glass  laser 
radiation.  During  the  course  of  study,  Q-switched  and  mode  locked  Nd:glass 
lasers  are  employed  as  excitation  sources  and  uniaxial  KDP  crystals 
of  various  crystallographic  orientations  served  as  nonlinear  optical 
media.  The  angular  dependence  of  relative  second  harmonic  intensities 
in  both  reflection  and  transmission  are  measured  and  compared  to  theory. 

The  theory  of  wave  propagation  in  a  medium  with  nonlinear 
susceptibility  and  criteria  of  conventional  phase  matching  and  noncollinear 
phase  matching  are  given  in  Chapter  II.  A  concise  account  of  pico¬ 
second  pulsewidth  measurement  using  phase  matched  second  harmonic 
generation  technique  is  also  given.  The  experimental  technique  including 
crystal  preparation  and  design  of  experiments  is  described  in  Chapter  III. 

In  Chapter  IV  detailed  experimental  results  are  given  and  discussed. 

The  proposal  and  brief  analysis  for  future  experiments  are  given  in 
Chapter  V. 


CHAPTER  II 
THEORY 


Introduction 


In  the  presence  of  very  high  intense  electromagnetic  radiation, 
the  medium  exhibits  an  electric  polarization  ^(r,t)  ,  which  is  in 
general  a  nonlinear  functions  of  electric  field  strength  E(r,  t).  The 
relationship  between  P(r,t)  and  the  applied  field  E(r,t)  can  be 


written  in  a  tensorial  form  as 


-+  m  •  -♦ 


P  =  Xi  E+  X2  :  E  E  +  x3  :  E  E  E  +  .  .  .  .  (2.1) 

where  x  1  s  are  nonlinear  susceptibility. 

If  the  applied  electric  field  is  harmonically  varying  in  time  t 

with  frequency  OJ  ,  then  the  second  term  on  the  right  hand  side  of  (2.1) 

is  the  nonlinear  source  term  at  the  second  harmonic  frequency 

*"  2  =  2  ««i  .  The  nonlinear  susceptibility  x  of  the  medium  will  give 

rise  to  a  polarization  at  the  harmonic  frequencies,  which  in  turn  radiates 

energy  at  these  frequencies.  The  effective  nonlinear  source  term  at 

the  second  harmonic  frequency  uu  =  2  0)  is  given  by 

2  1 

pNLS  =  pNLS(2(|i^  =  =  2  u’j)  :  E^E^  exp  i(k*" .  r-2  UJ^)  (2.2) 

T 

where  E  is  the  amplitude  of  transmitted  electric  field 
in  the  medium 

-♦ 

is- 

k  is  the  wave  vector  of  the  source  term  and 

Z’-Zk/ 

The  nonlinear  source  term  P  was  introduced  by  Armstrong 
.(U) 

et  al  .In  their  paper  the  nonlinear  susceptibility  tensor  y  was 
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shown  quantum  mechanically,  to  be  related  to  the  atomic  properties 
of  the  medium.  Furthermore,  they  also  showed  that  the  effective 
nonlinear  source  term  can  readily  be  incorporated  in  Maxwell'  s 
equations  for  the  nonlinear  medium.  Since,  during  the  course  of  study,  j 
large  number  of  harmonic  photons  are  created  by  the  interaction,  the 
problems  can  be  treated  classically  via  Maxwell'  s  equations.  In  this 
thesis,  we  shall  recast  the  theory  of  Bloembergen  in  a  form  that 
direct  comparison  with  the  present  experimental  data  can  be  readily 


made. 


A .  Wave  Propagation  in  a  Medium  with  Nonlinear  Susceptibility 


Maxwell'  s  equations  in  c.  g.  s.  units  for  nonlinear  medium  are 


V  •  D  =  4tt  n 


V  .  B  =  0 


vxS 

c  dt 


(2.3) 


^  x  H  =  —  D  +  —  J 

c  8t  c 


where  D  =  e  E  +  4tt  P 


NLS 


For  the  cubic  but  noncentrosymmetric,  nonconducting  crystal, 
the  wave  equation,  at  the  second  harmonic  frequency,  obtained  from 
(2.3),  is  given  as 


21  .  2-+NLS 


v  x  v  x  e  +  JL.e  2-  in 

2  2  7  2  2 

c  Bt  c  dt 


(20)  ) 


(2.4) 


2-+ 


2  j+NLS 


v2  £  . e(2w >  b__e_2  _  In  a  f 
2  2  2  "  2  2 
c  Bt  c  at 


(20)) 


Equation  (2.4)  is  similar  to  the  usual  linear  wave  equation  except 
it  is  augmented  by  a  source  term  on  the  right  hand  side.  The  general 
solution  of  (2.  4)  consists  of  the  solution  of  the  homogeneous  equation 
plus  the  particular  solution  of  the  inhomogeneous  equation.  The  general 
equation  for  the  transmitted  and  reflected  second  harmonic  field  created 
by  the  fundamental  field  incident  on  the  medium  from  the  vacuum  (see 
Fig.  2  -  t  )  are  given  as  , 

-T  A  T  -T  ->  4ttPNLS(4«,  2/c2) 

E2  =eTS2  expi(k2  .r-2Klt) - — - L_ -  * 

(k2  )  -  (k  ) 

-♦S..-+S  A 

FA  U  (k  .  p  “1  ~>s  -♦ 

P  -  - r.-1-  exp  i  (k  .  r  -  2u),t) 

L  no  1 


"*T  c  T  A  T  -*  'r  -» 

H2  =2^(k2  xeT)£2  exP  1  (k2  •  r  -  2tt)  jt ) 


(2.  5) 


-4TTPN^(4.i2/c2) 


(k2T)2  -  (kS)2 


c  "s  A  -*s 

2^. —  (k  x  p)  exp  i(k  .  r  -  2(1)^) 


p'  -  unit  vector  along  the  direction. 

The  reflected  harmonic  fields  are 


R  _  A  _  R  .  R  -» 
e2  ~  eR  S2  exp  i(k2  .  r  -  2-tJU  t ) 


(2.6) 


\  ’  R  A  .  R  R 

k2  X  eR  ?'2  exp  *  r  ”  2 


The  direction  of  the  wave  vectors  of  the  reflected  wave  k2  R  and 

*4  pT"’ 

the  homogeneous  transmitted  wave  k2  ,  as  well  as  the  polarization 

.  ^  R  T 

unit  vectors  e^nd  eR  and  the  magnitudes  ^  and  S2  have  to  be 

determined  from  the  boundary  conditions. 


i'l 


Z  z 


FIGURE  2  -  1 
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The  tangential  components  of  E  and  H  should  be  continuous 
everywhere  on  the  boundary  at  all  times.  This  requires  that  the  individual 
frequency  components  at  U  }  and  2^,  are  separately  continuous  across 
the  boundary  m  order  to  satisfy  this  condition  at  all  points  on  the  boundary 
simultaneously  .  Thus  for  the  fundamental  frequency 


k  1 =k  R -k  T 
lx  lx  "  lx 


(2.  7) 


The  polarization  of  electric  fields  at  2m  are  determined  by  the 


polarization  of  PNLS.  Also  it  is  required  that 


2k  T  =  k  S  =  k  R  -  k  T 
lx  2x  2x  2x 


(2.  8) 


Relations  (2.  7)  and  (2.  8),  in  short,  are  the  requirements  of  the 
conservation  of  tangential  component  of  momentum. 

L  -eneral  Laws  of  Reflection  and  Refraction.  Consider 


two  incident  plane  waves,  exp  i(k  1  .  r  -  »  t)  and  E  exp  i(k  \  r  - 


UJ2t)  ,  where  uj  2  4  U)  2  approaching  the  boundary  between  the  linear  and 
nonlinear  media,  from  the  side  of  the  linear  medium.  The  boundary  is 


defined  as  the  plane  z  -  0  as  depicted  in  Fig.  2-2.  The  angles  of 


incidence  of  the  two  ^V3.ves  3. re  P  ^  j  a  ^  rpi  i  .  , 

waves  are  ^  and  9 2  .  The  planes  of  incidence 


make  an  angle  cp  with  each  other, and  cp  varies  between  0  and  2tt  . 
Choose  the  z  and  y  direction  of  the  coordinates  system  such  that 


,  .  1  _  ,  i 

<-ly  2y  '  1x1  §enera1'  waves  at  all  sum  and  different  frequencies 


ml  %  -  m2  ^2  WiU  emanate  from  the  boundary  (n^  and  m2  are  integers). 
However,  in  case  of  second  harmonic  generation,  the  sum  frequency 


^  3  ^  +  w  2  2  w  2  ^  be  considered  explicitly. 
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A  necessary  and  sufficient  condition  for  the  requirement  that 
the  boundary  conditions  will  be  satisfied  simultaneously  at  all  points  in 
the  plane  z  =  0  is  that  the  x  and  y  components  of  momentum  wave 
vector  remain  conserved.  For  the  sum  wave,  this  gives  the  conditions. 


k  R=  k  T  =  k  S=  k  T+  k  T-  k  i  +  k  1 
3x  3x  3x  lx  2x  klx  +k  2x 


(2.9) 


R  .  T  .  S  .  T  .  .  T 


k3y  =k3y  =  k  3y  =  Sy  +k2y  =kly1+k2y1=0 


Equation  (2.9)  suggests  that  the  inhomogeneous  source  wave,  the 

homogeneous  transmitted  and  reflected  waves  and  the  boundary  normal 

all  lie  in  the  same  plane  which  is  xz  plane  called  "plane  of  sum 

reflection",  as  indicated  in  Fig.  2-2.  The  propagation  of  the  inhomogenous 

wave,  proportional  to  P  (0!  )  ,  is  given  by  exp  i{(k  L  +  k  ).r-(u)+uu  )t). 

J  12  12 

Its  angle  with  the  normal  into  the  nonlinear  medium  9^  is  determined  by 


sin9  S  =  |  k  T  +  k  T  1  '  '•  T 

3  lx  2x 


I  /  Ik/  +  k2T  | 


(2.10) 


T  T 

The  wave  vectors  k  ^  and  are  given  by  Snell*  s  law  for 

refraction  in  the  usual  linear  case.  From  Fig.  2-2  the  relationship 


between  k 


T  ;*  i 


’  ki  ’  k2  *s  S^ven  as 


|fc3T  l2  si„29  3T  =  |k  R|2sin29,R 


Ik^l  sin  9j'+  Ik  2  '  I  sin  9^  +  2  |kj  I  Ik^^.l  sin9|  sine^  cos 


CD 


Since 


e 


2 

c 

2  »  equation  (2.11)  can  be  written  as 

w 


(2.11) 


M  i„: 


. ..  u,  in  u. 


P 


I 

I 

I 

f 

I 

I 

I 

l 


e3Tc«32  sin293T=  c3Rt,32sin293R 


R  2  2  i  n  7  7  • 

el  s^n  +  62  ^2  s^n  ^2 

■  9  /  R  R.  1/2  .  i  i 

+  2  S1  e2  ^  lU)2  sin9i  sin^2  COS  ®P  (2.12) 

2  uu 2 

Since  k  =  e—  ,  and  also  from  (2.  9)  we  can  arrange  (2.  12)  as 

c 


rSsin2fl  S-  T-2qT  R  .  2  R 
e  Sin  9  -  e  sin  9  =  e  sin  9  „ 


(2.  13) 


Equation  (2. 13)  gives  directional  relationship  between  k  R  k  T 

-»  s  3  ’  3 

and  k  ^  .  This  relationship  is  the  nonlinear  counterpart  of  linear  Snells 

law. 

The  example  of  second  harmonic  generation  follows  from  equation 
(2. 12)  if  we  put  9^  =  9^  and  cp  =  0  ,  and  e.R  =  e  R 
This  will  lead  to  Snell1  s  law  for  nonlinear  case 


2  ,  and  U)  ^  =  uu  ^  =  uu  , 


T  .  2  t 
e3  Sln93 


R  .  2  i 
e1  sin  9 


or  in  general  (with  the  combination  of  (2.13)). 

/~R  .  „i  /  T  .  T  /  R  r 

VSj  sin  9  sin  9  =vS,  sin  9 


S  •  nS 

e  3  sin  9 


(2.14) 


I 

I 

I 

I 


The  relationship  (2.  12)  can  be  generalized  further  to  the  case  of 
two  beam  spatial  mixing  (TBSM)  second  harmonic  generation  where  two 
fundamental  beams  are  incident  from  the  opposite  side  of  the  face  normal 
(z  axis).  Here  we  have  9,1  ‘  1 


*1  ~  ^2  ’  ^  =  ^0  •  From  the  equation  (2.  12) 


and  (2.13)  we  get 


eT  =  es  =  9®  =  o  . 


This  means  that  all  thr 


(2.15) 


ee  second  harmonic  waves  travel  along 
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facing  normally  inward  for  and  kS  ,  and  outward  for  k*  with 
respect  to  the  boundary. 

From  the  above  analysis,  one  notices  that  the  conditions  of 
conservation  of  tangential  component  of  momentum,  equation  (2.9),  are 
very  general.  They  can  easily  be  used  to  derive  the  directional  relation¬ 
ship  for  higher  harmonics.  They  hold  regardless  of  whether  the 
harmonic  radiation  is  of  dipolar,  electric  or  magnetic,  or  quadrupole 
origin. 

2 •  Polarization  and  Intensities  of  the  Harmonic  Waves.  In 

this  section  particular  attention  will  be  devoted  for  the  polarization 

and  intensities  of  harmonic  waves  at  sum  frequency  uj  =  uj  +  m  =  2  m 

3  12 

The  nonlinear  polarization  induced  in  the  medium  is 


"*NLS  ^  -*  T  -*  T  r  T1  ^ 

P  («>z)  -  X  (2tt>)  :  Ey  E2  exp  ir(k1T  +  k_2T)  .  r-2u)tj  (2.16) 

The  refracted  waves  at  0)  1  and  U)  are  known  in  terms  of 
incident  waves  by  means  of  the  formulas  of  Snell  and  Fresnel  for  the 


incident  medium.  Therefore  the  nonlinear  source  term,  given  in 
Equation  (2.16),  at  U)  3  is  known.  The  angular  dependence  of 
itself  is  derived  from  the  transmitted  linear  waves  given  by  Fresnel 
equation.  One  must  take  proper  account  of  this  in  analyzing  the 
angular  dependence  of  harmonic  generation.  As  in  the  linear  case, 
the  second  harmonic  wave  at  U) ■  =  2  uu  with  electric  field  vector 
normal  to  the  plane  of  reflection  (E^  (2u>)),  as  defined  in  the  preceeding 
section^can  be  treated  independently  from  that  with  electric  field  vector 
in  the  plane  of  reflection  (E^(2U))) 
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2a.  Perpendicular  Polarization  (Ej_(2uu))  .  (E  =  E  , 
px  pz  ®).  Consider  the  case  where  the  nonlinear  polarization  is 
normal  to  the  plane  of  reflection,  i.  e.  PNLS  =  P  ^NLS  as  shown  in 
Fig.  2-3,  The  continuity  of  the  tangential  components  of  the  solution 


in  equation  (2. 5)  and  (2.6)  at  the  boundary  requires  in  this  case 

4ttP  NLS  ' 

E  -  FR  -  pT  x  1 

V  1  1  -  CrJ 


(2.17) 


Hx  =  _'/eRE-L  cos  9R  =  eTeI  c°s  9, 

.  . —  _  NLS 

4tt  /  e  P_l  cos 9  s 
+ — 2 - 

es  -  CT> 


(2.18) 


The  continuity  of  the  normal  components  of  D  and  B  follows 


automatically  from  the  equations  (2.17),  (2.18),  and  (2.12). 

R  T 

By  solving  for  E^  and  E^  from  equations  (2. 17)  and  (2.18),  one 

obtains  after  algebraic  manipulation. 

,NLS 


,R  '  4,tPL 
'1  ' 


ct‘es 


r  ^^TCOS,T-7es  cosqs  -I 

I  _  _ 

€t  cos  9t  +  ,/eR  cos9R  " 


(2.19) 


et  = 


.  t^nls 

4hP.  *»r,  -♦ 

J-  r  S  '  ~c 

-  I  exp  i(k  .r  -  2  out )  - - 


•/  ee  cos9g  +  7  eR  cos9r 


St"  €  S 


,/eT  cos9 T  +  y  eR  co s9r 


-+T  -♦  -| 

exp  i  (k  ,  r  -  2  U)t )  ; 


(2.  20) 


Since 


T  _  uu  r  / — 


A 


k  "  k  “  c  |  J  ®s  cos9s  "  V  CT  c°s9t  1  Z 


(2.21) 


Equation  (2.20)  can  be  transformed  into  a  single  plane  wave  with  a 


propagation  vector  k  ,  with  an  amplitude 


T  R  ^pNLS  ^ 

E  =  El  +  i~T~  fexp[  c(^~s  co^s-V^coseT)z]-l!  (2 . 22  a) 


For  the  small  values  of  z  which  satisfy  the  condition 


c  z'  '/eQ  cosi  -  </  S  tcosQ  '  «1 


(2 .  2  2  b  ) 


The  expansion  of  right  hand  side  of  equation  (2.22)  will  give  ET 
as  a  linear  function  of  z.  Then  the  intensity  grows  as  z2  from  the 
boundary.  This  is  precisely  the  effect  of  harmonic  generation  in  the 
volume  of  an  infinite  medium,  discussed  byArmstrong  et  al(11). 

2b.  Parallel  Polarization-  (E  =  P  NLS  =  0)  .  The  nonlinear 
"*NLS  ^  y 

polarization  P  has  only  x  and  y  components  as  indicated  in  Fig.  2-4. 

It  will  be  advantageous  to  describe  the  nonlinear  polarization  in 
the  plane  of  reflection  by  its  magnitude  PnNLS  and  the  angle  * 
between  its  direction  and  the  direction  of  propagation  in  the  source  k9*. 

The  continuity  of  the  tangential  componenets  at  z  =  0  requires  that 
from  equations  (2.  5)  and  (2.6). 

.  _  NLS  . 

r  T  4tt.F  sincycos9 

Ex  =  •  En  cos9r  =  V  co.eT  +— ue  - * 

es  eT 


,  _  NLS 

4rr  ll  coscy  sin  9^ 


(2.23) 


,  T3  NLS  . 

_  r  rp  *tt  i-1  sin** 

Hy  =  ”  ^  eRE  =--^Ten  -^S— ^ - 

eS  eT 


(2.24) 


By  solving  for  E  from  equation  (2.23)  and  (2.24),  one  obtained 


wave  E..  after  arranging  terms  as 


the  amplitude  of  the  reflected 


R 


A  NLS  .  .2 

4ttP.,  sin9  sm  9  sin  (w  +  9_  +  9  ) 
_ 5  I  T  S 


n 


11  eRSinqRSin(9T+0S)sin(q'r+9-)  cos(0t-9„) 


(2.25) 


T  R 


T  R 


With  the  introduction  of  the  angle  8  between  E  T  and  the  direction 
of  propagation  k  the  transverse  component  of  the  total  transmitted 
wave  will  be  given  with  the  help  of  equations  (2.25),  (2.26),  (2.23)and 
(2.24),  as 

T  4rrp  sin9  sin9  sin((y+9_+8_) 

E„  sinft  =  - - - - - I _ I _ 

11  eTsin(9T+9s)  sin(9T+9R)  cos(9T-6R) 

.  _  NLS  .  .  A 

4TT  11  Sln»  sin9gCos9g 

e  Tsin(9T+9s) 

+  4rrP11NL5coSasin(9T-9s)  expi  )} 


NLS  U) 

+  4nPu  sing/ cos(9T-9s)  exp  if-  z (/egCosBg-./  eTcos9T)}  -  1 


®S  "  ET 


(2. 26) 


The  longitudinal  component  of  the  electric  field  vector,  parallel 


to  ,  can  be  written  in  the  form 


-  T  Q  —  A  NLS  .  .  2 

11  COS^  ~  4ttP11  sln{y  Sln  9s-cos^cos(9T-9s)  sin(qT+8g) 

(JU 

exp  if-  z(,Aiscos9s-/lT  cos9t)] 

S T  sin(9T+  9^ )  (2.27) 

Nonlinear  Brewster's  Angie.  For  the  case  where  p^^P^  lies 
in  the  plane  of  reflection, the  electric  field  of  the  reflected  harmonic  wave 
is  given  by  equation  (2.25).  Equation  (2.25)  reveals  the  existence  of 
Nonlinear  Brewster' s  angle  for  harmonic  waves,  when  E^^O  .  Or  in 


LINEAR 


FIGURE  2  -  5 


other  words 


p 

EU  =  0  imPlies  +  eg  +  *  =  0  ,  tt  (2.28 

This  condition  implies  that  the  nonlinear  polarization  p^NLS  is 
parallel  to  the  direction  of  propagation  of  the  harmonic  wave  in  the 
nonlinear  medium  and  hence  it  can  not  radiate.  This  nonradiating 
wave  upon  refraction  back  into  the  linear  medium  would  otherwise  give 
rise  to  the  reflected  ray  in  the  direction  9R  .  This  is  shown  in  Fig.  2_5 
In  the  present  work,  several  KDP  crystals  of  rectangular 
parallelepiped  shape  are  used  as  the  nonlinear  media  from  which  second 
harmonic  waves  are  created  and  furthermore  the  optically  dense  fluid 

(1-Bromonaphthalene)  is  employed  as  the  linear  medium.  It  would  be 

advantageous  not  only  to  depict  wave  vectors  of  fundamental  and  second 

harmonic  waves  at  the  boundary  of  KDP  crystal  immersed  in  the 

optically  denser  liquid  (  Fig.  2.6  )  but  also  to  recast  those  fundamental 

equations  given  in  the  previous  section  into  appropriate  forms  such 

that  they  are  suitable  for  experimental  comparison. 

According  to  geometry  of  Fig.  2-6,  the  angles  0  ,9  and  9 

R  S  T 

of  the  reflected,  transmitted  source  and  transmitted  homogeneous  waves, 
respectively,  are  given,  with  the  a  id  of  equation  (2.  14)  and  n  =  , 

by 

niiq(u.  )  sine.  =  niiq(2'w)sin9R  =  n(a) )  sin9s  =  n(2uu)sin9T  (2.29) 

The  index  without  subscripts  refers  to  the  KDP  crystal. 

Since  the  liquid  is  optically  denser  than  the  KDP  crystal,  then 
there  exist  critical  angles  under  which  th 


e  fundamental  and  second 


harmonic  waves  are  totally  reflected.  Those  are  given  by 


sin9  (uu  )  =  nluu  )/n  (ur)  <  1 
cr  l.iq 

sin9  (2uj  )  =  n(2(U  )/n,  .  (<t) )  <  1 

cr  1  lq 
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(2.  30a) 
(2.  30b) 


The  components  of  the  harmonic  polarization  PNLS  along  the  cubic 
axes  of  the  nonlinear  KDP  crystal  are  given  in  terms  of  the  fundamental 
field  components  at  each  point  inside  the  crystal  by 


^  NLS  ,  Mr  r  'r 

Pz  (2U))=X36  \  (^)EyT(U) 


(2.  31a) 


p  NLS  NLS 

x  U,')  ’  Py  can  be  obtained  by  cyclic  permutation  of 

equation  (2.  31a ). 

Equation  (2.  31a)  can  be  expressed  in  terms  of  the  amplitude  E  of 

o 

the  incident  fundamental  wave  by 


^NLS  NL  L  2 

P  2'U)  =  X  3  6  ^FT  V 


(2.  31b) 


Where  q  is  a  geometrical  factor  which  depends  on  the  orientation 

of  the  fundamental  field  vector  and  nonlinear  polarization  component  with 

respect  to  the  crystallographic  cubic  axes  of  the  KDP.  The  linear 

Fresnel  factor  describes  the  change  in  amplitude  of  the  fundamenta: 

wave  on  transmission  at  the  crystal  surface.  For  the  laser  polarization 

perpendicular  to  the  plane  of  incidence,  it  is  given  by 
t  2  cos9 


cos9.  +sin9  (uu)  cos9 
i  cr  S 


(2.  32a) 


and  for  the  laser  polarization  in  the  plane  of  incidence  it  is 
T  2cos9. 


sinQ  (W)  cos9.  +  cos0  „ 
cr  l  S 


(2.  32b) 


The  nonlinear  polarization  PNLS  j 


is  the  source  of  the  three 


armonic  waves.  The  amplitudes  of  electric  field  of  second  harmonic 
waves  can  be  expressed  in  terms  of  PNLS  as  follows: 


E  (2oj)  =  4rr  PNLSp  NL 
K  R 

Es(2ou )  =  4tt  PNLSpgNL 

E  t  i) )  =  4tt  PNLSF  NL 
T  T 


(2.  33a) 


(2.  33b! 


(2.  33c) 


.  _NL, 

where  F  s  are  the  nonlinear  Fresnel1  s  factors. 

From  equations  (2.19),  (2.20),  (2.  22)  and  (2.33),  the  nonlinear 

Fresnel  factors  for  the  case  of  second  harmonic  polarization  perpendicular 
to  the  plane  of  reflection  are  given  by 


F  NL  e  —  1  .  At  r./yo.9T  -,/c;  coseg 

R.l  eT  “  ®s  */  Sg  cos9T+ye^cos0R 


)  <2' 


S.  1  “eT-ss 


(2.  34b) 


®T  "  eS 


(2.  34c 


For  the  case  of  second  harmonic  polarization  paralle  to  the  plane 
of  reflection  the  nonlinear  Fresnel  factor  can  be  obtained  from  equations 
(2.23),  (2.24),  (2.25)  and  (2. 33).  They  are  given  by 


NL  _  sinqgsin  9Tsin(^ +9T  +9g) 

R,U  eRSin9Rsin(9T+9R)cos(9T"9R)sin(9T+9s; 


NL  _  sing 


SS-eT 


j  5  sing 

®T  ®s"eT 


(2.  35a) 


(2.  35b 


(2.  35c) 


2  Q 

The  time  average  second  harmonic  power  carried  by  the  harmonic 
beams  is  given  by  the  real  part  of  Poynting  vector  multiplied  by  the 
cross-sectional  area  A  of  the  respective  beams.  The  intensities  of 
harmonic  beams  are. 

IR.S.TP“,  =  S-^.s.T  1  ER,S.Tt2",|2AR,S,T  (2'36> 


where  subscript  R,S,  and  T  refer  to  reflected,  inhomogenous  and  ho 


mo- 


geneous  harmonic  beams  respectively.  A  is  the  cross-section 

R >  o,  T 


area  of  the  beam  and  it  can  be  written  as 
.  d  d(  cos9 


R,S,  T 


R,S,  T 


COS0. 


(2.  37) 


where  d  d  '  are  rectangular  slit  which  defines  the  size  of  the  incident 
laser  beam. 

Since  the  intensities  of  harmonic  beams  cannot  be  observed  inside 
the  crystal,  it  is  preferable  to  relate  the  intensities  of  transmitted  beams 
after  they  have  reenterred  the  liquid  at  the  side  of  the  crystal.  The 
transmitted  power  from  the  crystal  to  liquid  is  given  by 


'R 


SS,T  1  S'T 


fLJ2 


(2.38a) 


where  f  is  the  linear  Fresnel  factor  for  transmission  at  the  interface 
S-  1 


for  crystal  with  right  angle  corner.  It  is  given  by 
,L  _  2 


l  +  f(e  -e_  _)  e(tt>)sin  2 9 .+ 1} 1 

R  S,  T  l 


(2.  38b) 


This  factor  is  always  close  to  unity,  and  f 


R 


=  1  by  definition. 


I 

I 

I 

I 

I 

I 

I 

1 
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I 
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I 

I 

( 
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Now  the  intensities  given  by  equation  (2.  36)  will  be  written  in  a 
rigorous  form  and  are  ready  for  direct  comparison  to  the  experiment. 


IR,S,T(2uj)=  (c/8tt)/T"r  |Eord  d/(4TTX36NL)2R2(f 

,2 


R,S,  T 


L  2 
)  x 


|fL|4  |fr,s,t^  cos9r,s,t<cos'1’)’1 


(2. 39) 


B •  Criteria  of  Optimum  Second  Harmonic  Generation 

In  this  section,  it  is  intended  to  point  out  general  criteria  for 
optimum  harmonic  generation,  in  particular,  second  harmonic  generation 
In  the  experiments,  which  verify  the  nonlinear  optical  laws  for  second 
harmonic  generation  (SHG),  it  is  anticipated  that  the  low  level  second 
harmonic  intensity  (SHI)  will  be  encountered.  Therefore  it  is  worth¬ 
while  to  utilize  the  condition  for  optimum  SHG  which  is  described  below. 

L  Phase  Matching  by  Birefringence.  When  a  fundamental 
wave  of  frequency  iu  propagating  in  a  direction  of  £  is  incident  on  a 
nonlinear  medium  of  piezoelectric  crystal,  the  nonlinear  polarization 

-*■  nls 

source  P  will  be  induced  in  the  medium  as 


PiNLS(2ur)=  E.(ur)  Ek(uu)  (2.40) 

The  nonlinear  polarization  PNLS(2u))  in  turn  will  radiate  SHI  in 
transmission  with  electric  field  E(2tu)  which  is  given  by 

2fjj 


(2.41) 


where 


A  k  =  k  r-jo)  -  2k  (w) 

2  ou .  2  uu  uj . 

=  —  (n  -n  ) 


(2.42) 


Consequently,  the  intensity  of  SHG  will  be  given  by 


A  I?  -> 


t/o  \  ,  ^0>  2  [  sin  2 

1(2  tw )  «  (x  )  ' 


.  r 


2 


1 


(2.  43; 


where  I  is  the  intensity  of  the  incident  laser, 
o 

One  possible  method  of  increasing  I(2(JU)  ,  described  in  equation 
(  2.  43  )  is  to  set  Ak  =  0  .  '  This  condition  can  be  achieved  by  the 

utilization  of  double  refraction  of  an  uniaxial  crystal,  e.  g.  KDP,  as 

/  q  q  v  /*}  A  \ 

pointed  out  independently  by  Giordmaine  and  Maker  et  al  .  This 


condition  is  called  phase  matching  which  is 


m  _  n  -  2t)U  /  2(U  . 

Ak  -  0  -  —  (n  -  n  ) 
c 


(2.44) 


Under  this  condition  both  fundamental  and  second  harmonic  beams 

travel  in  step  inside  the  crystal  with  the  same  phase  velocity.  However, 

when  Ak  4  0  ,  equation  (2.41)  shows  E(20)  )  to  vary  periodically  as  a 

function  of  distance  r  from  the  crystal  surface.  The  period  of 

variation  is  determined  by  the  term  so  called  "coherence  length"  i  , 

c 


which  is  given  by 
2ti 


2tt 


X 


c  Ak  ,  2W  w  20)  0) 

k  -2k  2  (n  -n  ) 


(2.45) 


The  physical  interpretation  of  f  is  that  it  is  the  maximum  crystal 

length  which  is  useful  in  producing  SHG.  It  is  noted  that  under  phase 

matching  condition  (Ak  =  0)  f  is  infinite,  In  practice,  however,  due 

to  double  refraction  phenomena  the  maximum  interaction  length  is 

limited  because  of  the  walk  off  effect  ^  ^  . 

The  diagram  showing  how  to  obtain  phase  matching  condition  in 

KDP  is  shown  in  Fig.  2-7.  In  this  diagram  the  angle  9  ,  with  respect 

m 


I 

i 


to  Z  (optic)  axis  is  called  the  phase  matching  angle.  For  KDP,  if  the 


fundamental  beam  is  launched  along  ^  as  an  ordinary  ray,  the  second 
harmonic  beam  will  be  generated  along  the  same  direction  as  an 
extraordinary  ray.  The  angle  9m  is  determined  by  the  intersection 
between  the  sphere  corresponding  to  the  normal  index  surface  of  the 
ordinary  beam  with  the  ellipsoidal  surface  corresponding  to  the  normal 
index  surface  of  the  extraordinary  ray.  Th 


.  2 
sm  9 


<n  V2-(„  2V 

_  o _ o 

m  ,  2io  v-2  ,  .-2 

(n  )  -  (n  ) 

e  o 


e  angle  9^  is  given  by 


(82) 


(2.  46) 


Taking  careful  consideration  to  equation  (2.  43)  the  criteria  for 
increasing  second  harmonic  generation  can  be  drawn  as  following. 

1.  By  using  laser  systems  that  provide  high  peak  power  or  large 

intensity  I  . 

o 

2.  By  selecting  proper  nonlinear  crystal  that  possess  large 
value  of  nonlinear  susceptibility  y  2a). 

3.  By  utilizing  phase  matching  condition  (Ak  =  0)  as  described 

above. 


2'  if.oncollinear  Phase  Matching.  According  to  the  previous 
section,  phase  matching  results  in  enhanced  optical  second  harmonic 
generation  (SHG).  In  certain  bifringent  crystal  such  as  KDP  this  is 
accomplished  by  making  the  birefringence  („*"  -  n^)  equal  to  the 
dispersion  (nQ  -  nQ  )  at  phase  matching  angle  9m  as  shown  in 
Fig.  2-7.  However,  SHG  may  still  be  limited  by  double  refraction(11’  14  ] 
which  is  describable  by  an  angle  0  between  the  poynting  vector  of 


^  r 


extraordinary  harmonic  wave  and  the  ordinary  fundamental  wave  -phe 

(35) 

angle  p  is  given  by 


where  9  is  the  angle  between  propagating  vector  of  homogeneous  second 
harmonic  beam  and  the  optic  axis. 


Consequently,  the  harmonic  wave  separates  from  the  fundamental 

for  a  beam  of  finite  diameter.  The  limitation  of  SHG  from  a  gas  laser 

(3  5) 

due  to  double  refraction  shows  the  reduced  .  production  of 


second  harmonic  beyond  the  aperture  length  £  which  is  given  by 


£  =  ,/tt  — 

a  0 


(2.48) 


where  w°  is  the  minimum  Gaussian  spot  size. 

Increased  SHG  will  result  from  £  m  when  p  =  0  .  However 

the  divergence  of  the  fundamental  beam  will  then  limit  SHG  instead  of 

the  aperture  length  £^  (or  double  refraction).  Phase  matching  with 

0=0  has  recently  been  achieved^42'  in  LiNbO  with  9  =  90°  by 

3  m  y 

varying  the  birefringence  and  dispersion  via  temperature  tuning. 

Another  phase  matching  method  with  p  -  0  is  via  noncollinear 
phase  matching.  This  geometry  is  handy  for  application  of  ultrashort 
pulse  measurement.  However,  one  should  point  out  that  even  when  p  is 
zero  as  in  this  case,  the  interaction  length  is  not  infinite.  One 
normally  does  not  adopt  this  scheme  for  maximum  harmonic  production. 

The  noncollinear  phase  matching  employed  in  the  present  work, 
requires  all  three  wave  vectors  k  1  s  forming  a  closed  loop.  The 


condition  is  given  by  several  authors^  '  ^  as  shown  in  Pig. 


(2uu )  -  k'j(ai)  +  k  '  (a) ) 


(2.49) 


where  kj  and  k  ^  are  wave  vectors  of  two  crossing  fundamental  beams, 

-4 

and  k  2  is  the  wave  vector  of  second  harmonic  beam.  According  to 
equation  (2.49),  this  vector  relation  can  be  satisfied  only  if 

"4  -4 

^2  ^  ^  ^  *  This  means  that  the  birefringence  must  exceed  dispersion. 

This  condition  is  applicable  to  KDP  crystal. 

(D 

n 

By  using  equation  (2.49)  and  relationships  k  ,  (uu )  =  -2—  •  ^ 

2u>  1  c 

i  ,  ne 

and  k.2(UJ)  =  c -  ‘  2  ^  °ne  Can  get  from  Fig-  Z~  the  conditions  for 

noncolinear  phase  matching  of  the  normal  component  : 

U).  2  u) 

%  n  o  n 

— —  .  U3  cos  (9  +9  )  + -  .uu  cos  (9  -9  )  =  — -  .  2tt)  cos  9  , 

c  13c  Z  5  c  i 


(2.  50a) 


For  the  tangential  component  : 


no  no  n 

— — -  •  tu  sin(9  +  6  )+ - u )  sin(9  -  9  )  = - UJsin0o  (2.50b) 

c  ioc  2  3c  3 

For  two  beam  spatial  mixing  (TBSM)  in  second  harmonic  generation, 

one  can  select  the  coordinates  system  and  crystal  orientation  such  that 

some  parameters  in  equation  (2. 50a)  and  (2.  50b)  will  be  eliminated. 

Such  configuration  is  illustrated  in  Fig.  2-9  .According  to  equation 

(2. 11)  and  the  symmetry  in  Fig.  2-9.  the  net  tangential  component  of  the 

■*  -♦ 

vector  kj  +  k  '  becomes  zero.  This  in  turn  makes  the  tangential 

-4 

component  of  k2(2oj)  zero,  i.  e.  9^  =  0  .  Under  this  condition  the 

-♦ 

propagating  vector  k  (2UJ)  will  be  in  the  direction  normal  to  the  optic 


(z  )  axis .  Equation  (2.50a)  ,  with  9^  =  0 

«J  U) 

n°  no 

oj  cos0  +  -  {ju  co s 0 _ 

c  I  c  2 


By  symmetry  of  Fig. 2-9,  9  =  0  =  cp 

1  2  m 

«>■  2uu 

n  cos  cp  =  n 
o  me 


,  can  be  reduced  to 

2  uu 

ne 

=  -  •  2oj  (2.  51) 

c  ' 

then  (2.  51a)  turns  out  to  be 

(2.  52) 


Equation  (2.52)  is  the  condition  for  noncollinear  phase  matching  for  the 


experimental  configuration  shown  in  Fig.  2-9.  Having  known  the  phase 

matching  angle  cp  inside  the  crystal,  the  phase  matching  angle  9* 

m  PM 

for  the  incident  beam  can  be  computed  via  Snell1  s  law  given  by  equation 
(2.14). 


9 


i 

PM" 


.  -lr  o  ~i 

sin  n  sin  Cp 

.  uj  mj 


(2.  53) 


(taking  n  .  =1.0), 

air 

C.  Theory  of  Picosecond  Pulsewidth  Measurement 

In  this  section,  the  theory  involving  the  utilization  of  nonlinear 

optical  methods  for  picosecond  pulsewidth  measurements  will  be  discussed.  Two 
photon  fluorescence  (TPF)  method  has  been  used  for  the  measurement  of 
picosecond  pulsewidth  by  several  authors (64’  ?2’  74\  However, 
particular  treatment  will  be  emphasized  on  harmonic  generation 
method,  which  is  utilized  in  this  experiment. 

'  nee  KDP  is  a  nonlinear  crystal  from  which  phase  matchable  SHG 
can  be  produced,  then  it  is  very  convenient  to  use  the  square  law 
intensity  characteristics  of  KDP  to  measure  the  second  order  auto¬ 
correlation  function  of  the  laser  light  which,  in  this  case,  is  a  picosecond 
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pulse  from  Nd:glass  laser.  In  the  previous  experiments  (62  ’  75,  78),  the 
harmonic  intensities  resulted  from  the  spatial  and  temporal  overlapping 
in  the  nonlinear  media  were  weak  such  that  they  could  be  detected  only 
Ly  means  of  photomultiplier. 

In  the  experiment  KDP  crystal  of  a  special  cut  such  that 
noncollmear  phase  matching  can  be  achieved,  is  served  as  the  source 
of  SHI.  Therefore,  it  is  anticipated  that  there  will  be  high  SHI  produced 
and  the  measurement  of  the  pulsewidth  can  be  done  by  means  of  both 
photomultiplier  and  photographic  methods.  The  theory  behind  this 
measurement  method  for  second  order  auto  correlation  is  given  below. 

l’  ?£cond  Order  Auto- Cor  relation.  The  SHI  radiating  from 
the  KDP  crystal  represents  the  variation  of  the  auto-correlation  function 
with  respect  to  the  degree  of  spatial  and  temporal  overlapping  of  the 
light  pulse  onto  itself.  This  can  be  performed  by  using  a  beam  splitter 
and  two  100<7,  reflectors  reflecting  the  two  beams  into  the  KDP  crystal 
where  in  spatial  and  temporal  overlapping  occurrs  .  The  schematic 
diagram  is  shown  in  Fig.  2-10.  The  incident  electric  field  is  E(r,  t). 

Then  the  resulting  electric  field  ER(r,t)  inside  the  KDP  crystal  is 
given  by 

■*  .-4  ■*  -4 

ER(r,t)  =  E(r ,  t )  +  E(r,t  +  t)  (2.  54) 

where  t  =  2nd/c  is  the  delay  time  of  one  of  the  reflected  beams,  and  2nd 
is  the  optical  path  difference  between  the  two  beams  .  The  output, 
resulting  from  the  overlapping  of  the  two  beam,  is  in  the  form  of 
average  second  harmonic  intensity  (SHE  in  the  direction  normal 


TO  PHOTOMULT 
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inward  to  the  crystal.  The  <'SHI\  can  be  expressed  in  terms  of  the  analytic 


(2.  55) 


signal  V  (t,r)  associated  to  the  real  field  E  (t,r)  as 

R 

(SHI)  ~  (  VR  (t,  r)  VR(t,  r)  VR''(t,r)  VRV(t,  r)\ 

'I' 

where  Vr  is  the  complex  conjugate  of  V  . 

The  analytic  signal  is  well  known  as  the  complex  representation 
of  an  oscillating  field.  The  real  part  of  V(t)  is  equal  to  the  light  field 
E(t)  and  the  imaginary  part  is  its  Hilbert  transform.  Since  in  the 
classical  limit  V  (t)  corresponds  to  the  quantum  mechanical 
anihilation  operator  a  and  V(t)  to  the  creation  operator  a^  .  One 
can  view  equation  (2.  55)  as  directly  related  to  the  quantum  mechanical 

probability  P  for  two  photon  absorption  in  the  two  photon  fluorescence  (TPF) 
method. 


P  ~  T/pa+  a+  aa  1  (2.  56) 

where  o  is  the  density  matrix  describing  the  light  field. 

By  using  equation  (2.  54)  and  (2.  55)  the  observed  6SHI}  at  a 

T 

certain  position  r  corresponding  to  the  delay  x  is  expressed  by  the 
incident  light  field  V(t)  as 

~  ((V- V  )  (V-V  )  (V _  v  '  )  (V  '  _  V* )  ) 

~  T  T  T  x 

*1,  j  , 

=  2(VVVV>  +  4(Y  V  V  '  V  ' 

T  T 

-  2(  (VV  +  V  V  '  )  (V  V  V  +  v'V  )) 

T  T  T  T 

-  (v  V  v  v  ""  +  v'1'  V  V "  V  )  (2  57) 

T  T  T  T  '  ’ 

where  V=  V(t)  and  Vt  =  V(t+  t)  . 

In  equation  {2.  57)  ,  the  first  two  terms  are  proportional  to  the 
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intensity  correlation  function  G(2),  known  as  Glauber'  s  second  ordei 
coherence  function'83 The  last  two  terms  describe  the  interference 
effects,  giving  spatial  variation  DroDorHnnal  ...  _  _  , 


i°n  Proportional  to  cos  w  j  and 

o 


cos  20) qt  ,  respectively,  where  0)  q  is  the  center  frequency  of  the 
light  field.  These  variation  are  averaged  out  to  be  aero  since  the 
average  value  of  coso,  qT  and  cos  2  V  over  and  interval  At 
corresponding  to  several  wavelength  vanish.  The  normalised  intensity 
'SHI1  now  expressed  by  intensity  correlation  function  C<2),t>  of 


the  incident  field 


(SHI)  =  (L9 _ (P)  +  4  C  \t  ) 


G(2)(0) 


(2.  58) 


where  G(  ' (r  )  is  defined 


G(2)(t)  =  V'"(t)  V'"ft+r)  ) 

(  v  (t)>  /V(t+T)  V  (t+  T)  ) 


(2.  59) 


The  bracket  indicates  the  time  average  or  ensemble  average.  If  only 

the  time  average  is  considered  .  we  can  use  y t,  .  VR,t  I  v*ft  I  and 
arrange  equation  (2.  59)  into  a  new  form  as 

I(t)  I(t+  r)dt 

? 


2  d 
r~c  ,  ,  -oo 

G  (T)  =  - - 


(2.  60) 


where  I(t)  is  actual  intensity  and  l  is  the  mean  intensity. 

In  the  experiment  the  KDP  crystal  is  oriented  in  such  a  way 
that  the  noncollinear  phase  matching  for  SHI  will  be  achieved  if  and 
only  if  the  two  beams  have  temporal  and  spatial  overlapping  inside  the 
crystal.  The  nonlinear  pola  risation  PNL?  will  be  c  rea  ted  and  then  it 
radiates  at  second  harmonic  frequency.  In  terms  of  those  applied 
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electric  fields  the  P  is  given  by 


T-.NLS 

P  ~  E(t)  .  E(t+  t  ) 


,NLS 


(2.61) 


Further,  P  can  be  described  by  analytic  signal  V  .  Using 
the  analytic  signal  V(t)  associated  with  the  real  field  E(t)  ,  the 
analytic  signal  V  is  directly  obtained  from 

o  H 


VSH  ~  V(t)  .  V (t+  T) 


(2.  62) 


By  using  I(t)  =  V(t)  v  '  (t)  and  equation  (2.  62 ).  The  average 
second  harmonic  intensity  f  I  >  produced  by  two  parts  of  the  same 
signal  V(t)  with  spatial  overlapped  and  time  delay  r  is  given  by 

^SH>  ~  V(t)  V(t+T)  V*(t)  V*(t+T)  ~G(2)(t  )  (2.63) 


In  this  case  the  second  harmonic  intensity  is  a  direct  measurement 
of  the  second  order  auto-correlation  function  G^(t)  .  Since  neither  one 
the  fundamental  beams  can  produce  SHI  in  the  normal  direction,  then 
there  will  be  a  background  free  auto  correlation  function  G^(t)  and 

there  is  no  contrast  ratio,  as  in  TPF  method,  to  be  considered. 

(2 ) 

Since  the  G  (t )  obtained  in  the  experiment  is  background  free 
function,  then  by  using  the  photographic  method  the  exposure  of  the 
film  will  reveal  the  pulsewidth  Qf  the  G(2)(  T  )  function  which  in  turns  gives 
the  direct  measurement  of  the  picosecond  pulsewidth.  Furthermore, 
the  two  fundamental  beams  are  obliquely  incident  on  the  crystal  as 
shown  in  the  Fig.  2-11,  Due  to  the  finite  size  of  the  two  beams,  there 
will  be  a  built  -in  time  delay  along  different  points  on  the  surface  of  the 
very  thin  crystal.  As  shown  in  Fig.  2-11,  at  the  central  point  of  the 
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overlapping  region,  the  light  pulse  of  two  beams  arrive  at  the  same 
time,  thus  there  is  no  time  delay.  At  the  point  x  away  from  the  center 
the  light  pulse  from  one  beam  will  arrive  at  time  At  earlier  then  the 
other  beam.  The  delay  is  given  by 

O  . 

At  _  2x  cos  (9O-01) 
c  /n 

.  2x  sin  9* 


where  c  is  the  velocity  of  light, 

9  the  angle  of  incidence 

n  the  refractive  index  of  the  linear  medium. 

For  a  beam  of  width  d,  the  maximum  delay  T  across  the  entire 


overlapping  region  is 


At )  =  —  sin  9* 

max  c 


where  n  =  1.  0  for  air. 


For  d  =  5  mm,  91  =  15°  the 


(2.  64) 


maximum  delay  At)  is  calculated 

max 


to  be  about  5  picoseconds.  This  total  built-in  delay  will  be  enought  to 
cover  a  picosecond  pulsewidth  which  has  a  theoretical  value  about  0.  33 


picosecond. 


In  order  to  cover  the  entire  picosecond  pulsewidth  measured  at  the 
base,  a  crystal  of  larger  birefringence  is  required. 

Such  a  crystal  will  provide  a  relatively  large  noncollinear  phase 
matching  angle  which  in  turn  gives  at  higher  value  of  A  T  )  in  (2.  64) 


Another  factor  that  will 


increase  At)^^  is  the  index  of  refraction  n 


of  the  linear  medium.  One  can  use  a  linear  optically  denser  liquid  called 
1-  Bromonaph  thalene  instead  of  air  in  the  experiment. 


CHAPTER  III 

EXPERIMENTAL  TECHNIQUE 


Since  the  present  experimental  work,  by  nature,  involves 
nonlinear  optical  interaction  in  KDP  crystal,  it  is  appropriate  to  concisely 
describe  the  experimental  arrangement s  and  techniques  including  the  laser 
systems 'hatwere  used  The  source  of  excitation  was  a  Nd:glass  laser 
which  was  operated  in  mode  locked  and  Q-switched  fashions.  During 
the  course  of  the  experiment  very  high  and  low  second  harmonic 
intensities  (SHI)  were  encountered  and  properly  detected  and 
subsequently  they  were  compared  to  the  theory.  Furthermore,  the 
experimental  results  concerning  polarization  properties  and  phase 
matched  condition  of  KDP  crystal  were  utilized  for  measurement  of 
picosecond  pulsewidth  of  Nd:glass  laser.  In  this  chapter  major 
experimental  preparations  and  techniques  involved  in  setting  up  and 
detection  second  harmonic  signal  will  be  given.  Since  a  Ndrglass  laser 
was  employed  as  the  excitation  source.and  also  to  facilitate  under¬ 
standing  the  picosecond  pulsewidth  measurement,  the  Nd:glass  laser 

operated  in  mode  locked  and  Q-switched  fashions  will  be  briefly  discussed 
below. 


A-  The  Neodymium  Glass  La 


s  er 


The  laser  rod  is  made  of  glass  in  which  the  Nd  ions  reside 


and  act  as  impurities.  The  energy  levels  involved  in  the  laser  trans- 


+  3 


iHon  the  edited  states  ofNd"  ion.  The  laser  acti 


ion 


take  place  from  the  upper  level  to  the  lower  level  4  ^ 
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4p 


is  approximately  1950  cm  above  the  ground  state  as  shown  in  Fig.  3-1. 
The  laser  is  classified  as  4  level  laser  system  from  which  the  emission 
wavelength  =  1.06pm  is  radiated.  The  fluorescent  linewidth  can  be 
measured  directly  and  ranges  around  300cm"1  (82  \  This  width  is 
rather  broad.  This  is  due  to  the  amorphous  structure  of  glass,  which 
causes  different  Nd+3  ions  to  '-see"  slightly  different  surrounding. 

This  causes  their  energy  splitting  to  vary  slightly.  Different  ions 
consequently  radiate  at  slightly  different,  frequenciespausing  inhomogeneous 
broadening  of  the  spontaneous  emission  spectrum.  This  larger  band¬ 
width  is  advantageously  utilized  for  riode  lock.ng  the  laser  since  it  can 

support  many  axial  modes  ,  N  ,  in  the  oscillation  .  For 

4 

Ndtglass  laser  N  is  about  10  . 

The  mode  locked  Nd:glass  laser.  The  mode  locked 

Ndtglass  laser  used  in  the  experiment  was  a  Korad  Kl-sysiem  as  shown 

in  Fig.  3-2.  It  is  consisted  of  Ndrglass  rod  having  a  Brews'er  -  Brewster 

configuration  to  avoid  reflection  losses.  The  laser  rod  is  of  diameter 

1.  25  cm  and  of  length  20cm.  The  two  coated  dielectric  mirrors  have  at 

1.  06  u  m  reflectivities  of  100%  and  65°i  .respectively,  and  form  a 

laser  cavity  having  75  cm  length.  In  addition  the  contact  dye  cell  of 

0.  078  cm  thick,  is  attached  to  the  mirror  M^lOO^,  reflectivity).  The 

cell  contains  Kodak  9860  dye  solution  dissolved  in  dichloroethane  .  The 

solution  serves  as  passive  nonlinear  absorber  which  is  required  to 

achieve  mode  locked  pulse.  Since  the  inhomogeneous  linewidth  (  Avlofthe 
12 

laser  is  3x10  Hz  ,  then  it  is  expected  theoretically  that  the  laser 


pulse  will  be  of  0.  33  picoseconds  duration. 
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To  generate  ultra  short  pulses  from  the  system,  all  the  laser 
modes  (axial  modes)  falling  within  the  laser  linewidth  must  be  coupled 
together.  This  can  be  achieved  by  introducing  the  cell  containing  passive 
nonlinear  absorber  as  described  above.  Despite  of  its  wide  use,  passive 
mode  locking  is  not  well  understood.  However  the  mode  locking  by 

using  saturable  absorber  can  be  explained  on  a  qualitative  basis  in  the 
time  domain  as  follows: 

The  laser  can  be  viewed  as  a  quantum  mechanical  oscillator 
building  up  from  spontaneous  emission  noises  of  various  amplitude  and 
intensity.  The  saturable  absorber  (dye  solution)  has  a  nonlinear 
abso  rption  characteristic  which  for  low  tight  intensities  acts  as  a  strong 
absorber  and  for  high  light  intensities  is  transparent.  When  the  laser 
starts  to  build  up,  the  low  amplitude  portions  of  the  amplified  fluctuating 
spontaneous  emission  noises  are  discriminated  against  the  higher 
amplitude  portion  because  of  the  nonlinear  absorber.  As  a  consequence, 
the  lower  intensity  portion  of  the  pulse  is  cut  off  and  the  peak  portion  is 
allowed  to  pass  through  the  absorber.  By  the  time  the  pulse  is  reflected 
back  from  the  mirror,  the  absorber  relaxes  to  the  ground  state  and 
once  again  the  pulse  is  sharpened  and  amplified  by  the  laser  rod.  After 
many  passes  thmugh  the  cavity  ,  the  pulse  is  narrowed  to  its 
limit,  which  is  the  inverse  of  the  oscillating  bandwidth  of  the  laser  system. 
The  output  will  be  a  train  of  picosecond  pulses  contained  in  a  single 
envelope.  The  characteristic  length  of  the  pulse  train  is  about  300-600 


nanoseconds. 


J 


The  mode  locked  laser  has  two  major  advantages  which  are 
employed  in  the  experiments.  First,  not  like  other  lasers,  it  has  a 
definite  phase  relationship  among  participating  modes  (axial  modes). 

This  causes  the  intensity  of  the  laser  output  to  have  a  temporal  coherence 
and  overall  less  fluctuation'87 >.  It  is  very  important  to  the  SHG  experiment 
aince  it  is  a  nonlinear  optical  process.  Inde<,d  ,  experimental  data  points 


are  found  to  exhibit  less  fluctuate 


Secondly,  the 


laser  system  of  this  type  has  very  high  peak  power.  Therefore  it  is  very 
useful  for  SHG  when  this  effect  is  extremely  small  as  in  certain  situation, 
e.  g.  Nonlinear  Brewster' s  angle  condition.Due  to  the  ultrashort  pulse 
duration  the  damage  threshold  to  the  crystal  will  be  high  ;  consequently 
crystal  will  not  be  easily  damaged  by  the  experiment. 

2>  Hi?  Q  switched  Ndtglass  Laser.  In  the  experimental 
investigation  for  two  beam  spatial  mixing  (TBSM)  and  for  noncollinear 
phase  matching  in  KDP,  the  Nd:glass  Q-switched  laser  was  utilized. 

Q -switched  pulses  were  obtained  by  means  of  Korad  Pockels  cell  model 
K-QS2  and  polarizer  stack  consi  sting  of  glass  plates  oriented  at  Brewster' s 
angle.  The  Q-switched  Nd:glass  laser  system  is  depicted  in  Fig.  3-3. 

The  Pockels  cell  and  polariz er  s tackwere  oriented  in  such  a  way  that 

their  polarization  axes  were  orthogonal  to  each  other,  A 

Pockels  cell  operation  can  be  briefly  described  as  follows  .  ft  serves 
as  an  optical  shutter  for  Q-switched  operation.  While  the  laser  rod 
is  being  pumped  and  the  laser  is  building  up  the  Pockels  cell  will  be 
passive  in  the  sense  that  it  blocks  or  prevents  the  laser  light  from 
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reaching  the  cavity  mirror  .  This  will  introduce  high  losses  into  the 
system  such  that  the  laser  oscillation  is  prevented  at  the  low  population 
inversions.  When  the  population  inversion  in  the  laser  rod  reaches 
maximum  value,  the  Pockels  cell  is  triggered  by  means  of  the  delay 
electronic  circuit.  At  this  moment  a  high  voltage  is  applied  across 
the  crystal  in  the  cell  such  that  via  electrooptical  effects  the  axis  of  the 
polarization  of  the  crystal  is  rotated  by  it  /2  .  The  laser  light  will 

passthrough  the  cell  .reaching  the  mirror  Mj  and  be  reflected  back.  Upon 

returning  from  mirror  Mj  and  passing  through  the  cell  again,  the 
polarization  of  the  light  is  suffered  another  tt/2  rotation.  The  total 
change  of  the  polarization  of  the  light  is  tt  .  Then  the  light  can  pass 
through  the  laser  rod,  the  polarizer  stag  and  reaches  the  mirror  M  . 

The  essential  role  of  polarizer  stag  is  to  enforce  the  polarization  of 
the  light  and  helps  the  Pockels  to  work  efficiently.  This  will  result  in 
raising  the  Q  of  the  system  to  a  very  high  value  in  a  short  time.  As  a 
consequence  the  system  will  oscillate  at  the  highest  population  inversion 

level.  Then  the  laser  will  lase  out  from  the  mirror  M.,  as  a  giant  pulse 

having  pulse  duration  about  30  nanoseconds.  The  Pockels  cell  for 
Q-switching  is  very  effective  and  the  Q-switched  pulse  is  reproducible. 
Besides  it  is  much  easier  to  operate  than  the  clumsy  passive 

dye  Q-switched  cell  used  in  the  early  day.  The  advantage  of  using  a 
Q-switched  pulse  for  two  beam  spatial  mixing  (TBSM )  experiment  is  that 
its  pulse  duration  is  so  long  that  the  temporal  overlapping  is  always 


achieved. 
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B.  KDP  Crystal 

The  nonlinear  crystals  used  in  the  experimental  investigation 
are  all  KDP  (Potassium  di  Hydrogen  Phosphate)  of  several  crystallographic 
cuts.  The  orientation  of  crystallographic  axes  of  the  crystals  are  given 
along  with  the  experimental  results  in  chapter  IV.  The  crystal  KDP  is  a 
tetragonal  crystal  belonging  to  the  symmetry  class  4  2m  .  It  is  piezo¬ 
electric  and  the  second  harmonic  polarization  in  it  is  given  by  the 

equations  (2.40)  in  Chapter  II  and  its  nonlinear  susceptibilities  has  been 

,  u  '84) 
given  elsewhere 

The  crystal  is  eminently  suitable  for  the  present  experiment 
since  it  is  phase  matchable  and  has  relative  high  value  of  the  nonlinear 
susceptibility.  Furthermore  it  ^s  transparent  at  the  fundamental  and 
second  harmonic  wavelengths  respectively.  This  will  enable  an 
investigation  in  transmission.  In  addition  to  possessing  the  intrinsic 
properties  described  above, its  linear  optical  properties  are  isotropic 


so  that  the  equations  developed  in  chapter  II  are  applicable  .  Because 

its  linear  index  of  refraction  is  relatively  low  ,  total  reflection 

from  it  is  possible  via  the  optically  denser  linear  fluid  1-  Bromo- 

/  o  2) 

naphthalene.  From  the  table  s  for  refractive  dispersion  its  indicies 
of  refraction  for  the  ordinary  and  the  extraordinary  rays  at  the  fundamental 
and  second  harmonic  wavelengths  are  deduced  and  given  as 


uu 

n  =  1. 4943 
o 

2  U) 

n  =1.  5131 
o 


2u> 

n  =  1.4708 
e 
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Where  (0  corresponds  to  the  fundamental  wavelength  of  1.  06pm. 

The  KDP  crystals  were  ordered  and  prepared  by  Gould  Inc. 
Tne  typical  dimensions  of  the  crystals  are  25x15x8  mm3.  The 
entrance  and  exist  surfaces  are  polished  optically  flat  to  \/5  at  the 
D-lme  of  sodium  light.  The  other  surfaces  are  see-through  polished. 
The  parallelism  between  opposite  faces  is  better  than  30  seconds,  and 
none  of  surfaces  are  coated. 


c-  Optically  Dense  Fluid 

In  the  experimental  investigation  of  a  phenomenon  so-called  total 
reflection  and  Nonlinear  Brewster's  angle  are  involved. The  investigation 
under  these  conditions  can  be  achieved  only  if  the  linear  medium  from 
which  the  laser  beam  is  incident  has  a  higher  index  of  refraction  than 
KDP.  To  observe  the  Nonlinear  Brewster's  angle  of  KDP  in  air  would 
require  an  angle  of  incidence  from  air  to  KDP  greater  than  90°.  This 
is  impractical.  Therefore.it  is  essentia!  that  the  linear  medium  in  contact 
with  KDP.must  possess  ahigher  index  of  refraction.  By  using  optically  denser 
fluid  1-Bromonaphthalene,  the  angle  of  incidence  which  corresponds  to 
the  Nonlinear  Brewster'  s  angle  is  in  the  vicinity  of  43.  0°  . 

In  the  experiment,  the  nonlinear  crystal  KDP  was  immersed 
in  the  optically  denser  fluid  1-Bromonaphthalene  which  has  larger  indices 
of  refraction  than  KDP  at  both  ®  andZui  The  fluid  is  transparent 
from  the  range  0.4  -  1.06  pm.  The  index  of  1- Bromonaphthalene  has 
been  tabulated  at  five  different  wavelengths  <85 '.  The  values  at  wave¬ 
lengths  of  interest  may  be  interpolated  by  means  of  the  Cauchy  relation186’. 
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This  gives  for  the  index  of  refraction  of  the  fluid 

n.  (uo )  =  1.  6260  n  (2(d’)  =  1.  6701 

l  1C1  liq 

F rom  Snell'  s  law  given  in  Chapter  II  one  finds  the  critical  angles 
for  total  reflection  of  the  fundamental  and  second  harmonic  beams  to  be 
cr^)  =  *^'78  c?cr"(2uJ  )  =  64.  76  ,  respectively. 

D.  Experimental  Arrangement 

In  the  present  study,  the  experimental  arrangement  can  be 
divided  into  two  major  parts.  One  involves  detection  of  SHG  in  the 
reflection  and  transmission  using  a  single  incident  laser  pulse  and  the 
other  is  two  beam  spatial  mixing  (TBSM)  setup  and  measurement  of  a 

picosecond  pulse.  The  experimental  arrangements  are  described  in  the 
following. 

Experimental  Arrangement  for  SHG  in  Reflection  and 

Transmission.  The  excitation  source  for  this  case  was  a  Ndrglass  laser 

operated  in  mode  locked  fashion.  The  laser  system 

cavity  mirrors  M  and  M  ,  a  contact  dye  cell  with  dye  solution  Kodak 

2 

9860  dissolved  in  dichloroethane  and  a  water  cooled  laser  head.  The 

+  3 

laser  head  consists  of  a  Nd  doped  glass  rod.  The  rod  was 

surrounded  by  a  helical  xenon  flash  lamp.  The  laser  system  is  shown  in 
Fig.  3-2  in  the  previous  section.  All  optical  alignment  was  performed 
with  the  help  of  He-Ne  gas  laser.  The  cavity  length  which  was  the 
optical  path  length  between  the  mirrors  and  M 2  is  about  75  cm. 

The  entire  laser  system  was  covered  with  a  wooden  box,  to 


TO 
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prevent  the  flash  lamp  radiation  during  the  pumping  period  from 

interfering  with  the  SHI  detection.  Furthermore  a  Corning  filter.no 

2-64  (C.  F.  2-64)  was  used  at  the  exist  of  the  wooden  box  to  filter  out  the 

residual  flash  lamp  light  in  the  lasing  direction.  The  experimental 

layout  which  will  be  described  hereafter  is  depicted  in  Fig.  3-4.  The 

fundamental  beam  was  then  partially  deflected  by  an  ordinary  glass  slide 

plate  into  a  ITT  photodiode  F-400  (S-l)  with  ground  glass  and 

neutral  density  filter  in  front  of  it.  The  pulse  from  the  ITT  photodiode, 

displayed  on  a  Tektronic  519  oscilloscope  had 

an  overall  rise  time  of  0.7  nanosecond.  In  turn  the  Tektronic  519 

oscilloscope  also  provided,  upon  registering  the  fundamental  beam,  the 

triggering  signal  to  the  Tektronic  557  dual  beam  oscilloscope,  used  for 

recording  SHI  from  the  monitor  channel  and  from  the  KDP.  A  piece  of 

z-cut  quartz  platelet  was  used  in  the  monitor  arm.  The 

laser  pulse  causes  the  quartz  to  generate  SHI  which  is  used  as  a 

reference  for  reducing  of  statistical  fluctuations  in  the  SHG  exoeriment 

The  use  ofa  monitor  channel  is  a  standard  technique  in  nonlinear 
(27,  28) 

experimentation 

The  reason  for  using  a  z-cut  quartz  platelet  quartz  in  the  monitoring 

arm,  rather  than  KDP,  was  that  its  surface  did  not  deteriorate  with 

time.  The  z-cut  was  chosen  mainly  because  of  the  convenience  in  not 

having  to  worry  about  the  orientation  of  the  laser  polarization  with 

(2  7 ) 

respect  to  the  crystalline  axis  .  If  the  z  axis  is  normal  to  the 
platelet,  the  second  harmonic  intensity  from  the  quartz  is  pro po rtiona  1  to 


5°> 
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Here  ly® )  and  yw  )  are  the  projection  of  the  laser  field  on 
the  X  and  Y  axes  of  the  quarto  crystal.  The  s econd  ha rmonic  intensity 
was  not  chanced  when  the  quartz  is  rotated  about  the  z  axis.  Another 
advantase  of  using  a  z-cut  was  that  the  crystal  was  far  from  the  phi 
matching  condition.  Thus  accidental  misalif 


iase- 


system  would  not  effect  the  second  ha 


ignment  of  the  monitoring 
rmonic  conversion  by  a  noticeable 


amount. 


CorninB  filter  CF  2-64  and  7-57  were  used  in  front  the  quartz 
platelet  in  order  to  ensure  that  only  the  fundamental  beam  could  strike 
the  quartz.  To  discriminate  the  fundamental  beam  from  the  admixture 
of  second  harmonic  signals  from  the  quartz  platelet,  a  copper  sulphate  soli 
lion  (CuSCM  was  employed  to  filter  out  the  fundamental.  The 
remaining  signal  is  then  transversed  through  a  Viard  Atomic  interferenc, 
filter,  having  transmission  peak  at  wavelength  5300A  corresponding  to 
second  harmonic  of  the  fundamental.  The  signal  was  collected 

ten  stage  Amperex  photomuItiplierfmodel  56  AVP.  The 
signal  from  the  photo  multiplier  in  the  monitor  channel  „as  fed  to  the 
upper  beam  terminal  of  the  Tcktronic  551  dual  beam  oscilloscope. 

The  mam  fundamental  beam  was  polarized  in  the  vertical 
direction  corresponding  to  the  ^0]  direction  of  tho  KDP  target. 

The  linear  polarization  of  the  fundamental  beam  was  achieved 
by  means  of  half  waveplate  and  Glan-Kappa  prism  which  were  coated  for 
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anti- reflection  at  the  wav-’  ngth  of  1.06pm  .  Before  the  fundamental  beam 
entered  a  liquid  cell,  it  transversed  through  Corning  filters  no.  2-64  and 
no.  7-56.  This  was  to  ensure  that  all  spurious  signals  generated  by  the 
laser  beam  from  various  optical  components  were  suppressed  and  only  the 
fundamental  beam  of  the  wavelength  )  =  1.  06  pm  would  be  incident  on 
the  KDP  crystal. 

The  KDP  crystal  was  mounted  onan  aluminum  target  holder 

which  was  connec  ed  to  the  angular  rotator  mounted  on  the  platform  above 

the  liquid  cell.  The  angular  rotator  was  Kinematic  model  RT  200,  with 

vernier  scale  permitting  variation  of  rotation  within  accuracy  of  0.  01°. 

The  fundamental  beam  before  entering  the  liquid  cell  was  regulated  by  a 

rectangular  slit  1  mm  wide  and  5  mm  high.  The  liquid  cell  had  a 

hexagonal  shape  with  a  circular  flat  fused  quartz  window  on  each  side. 

Each  window  subtends  an  angle  of  20°  .  The  cell  contained  optically 

denser  fluid  1  - Rromonaphtha lene  in  which  the  KDP  crystal  was  immersed. 

The  KDP  crystal  was  oriented  in  such  a  way  that  its  nonlinear  polarization 
^NLS 

source  P  ,  m  the  direction  of  optic  fz)  axis,  lay  in  the  plane  of 

reflection.  The  fundamental  beam  was  polarized  normal  to  the  plane  of 

reflection.  In  other  words, the  polarization  of  the  fundamental  beam 

was  m  the  direction  HTo]  with  respect  to  the  crystallographic  axes  of 
the  KDP. 


The  detecting  system  for  SHI  produced  from  KDP  was  mounted 
on  an  aluminum  arm  pivoted  underneath  the  liquid  cell.  The  axis  of 
rotation  of  the  arm  was  common  to  the  line  passing  through  the  center 
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of  the  cell  and  tangent  to  an  entrance  surface  of  the  KDP.  This  will 
ensure  that  when  the  arm  was  rotated  to  change  the  angle  of  incidence 

91  ’  the  SHI  would  bc  Properly  collected.  Furthermore  it  helped  to 
verify  the  directional  property  of  SHI  by  rotating  arm  slightly  off  to  both 
sides  of  the  expected  direction  of  the  SH  signal. 

The  detection  system  for  SHI  from  the  KDP  wa  ,  consisted  of  the 

following.  A  slit  4mm  wide  and  10  mm  high  was  used  to  separate  the 

reflected  or  transmitted  SHI  from  the  reflected  or  transmitted  fundamental 

beam  respectively.  The  slit  was  placed  about  50  cm  away  from  the 

liquid  cell.  For  the  case  of  total  transmitted  SHI  where  there  were  two 

transmitted  beams  of  homogenous  and  inhomogeneous  SHI,  the  larger 

slit  was  placed  at  15  cm  away  from  the  cell  to  allow  the  two  harmonic 

beams  to  pass  through.  Behind  the  slit  a  biconvex  lense  of  30  focal 

length  was  employed  to  ensure  that  all  the  harmonic  signals  would  be 

properly  collected  by  the  photomultiplier  which  was  about  75  cm  away  from 

the  liquid  cell.  Between  the  lense  and  the  photomultiplier  a  copper 

sulphate  cell  and  Baird  Atomic  interference  filter  with  transmission  peak 
o 

at  5300A  were  employed  in  the  manner  as  described  previously  for  the 
monitor  arm  channel.  Between  the  copper  sulphate  cell  and  the  inter¬ 
ference  filter,  a  sheet  of  polariod  was  used  for  verifying  the  polarization 
orientation  of  second  harmonic  signal  before  passing  through  the  inter¬ 
ference  filter.  The  photomultiplier  in  this  channel  was  identical  to  that 
of  the  monitor  channel  so  that  the  rise  time  of  the  two  channels  will  be 
essentially  the  same  and  time  correlation  for  SHG  could  be  easily  verified. 
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The  signals  from  both  channels  were  displayed  on  the  dual 
beam  oscillscope  and  photographed.  The  oscilloscope  was  triggered 
externally  by  a  triggering  signal  from  the  Tektronic  519  oscilloscope  which 
simultaneously  monitored.  To  generate  a  data  point,  we  took  the  ratio 
of  the  pulse  heights  corresponding  to  the  signal  and  monitor  intensity 
respectively  and  then  averaged  this  ratio  over  five  to  ten  laser  firings. 

This  was  essential  in  order  to  improve  the  quantum  statistical  fluctuations 
of  the  SHG.  The  data  points  were  plotted  and  compared  to  the  computed 
theoretical  curves. 

2.  Experimental  Arrangement  for  Two  Beam  Spatial  Mixing 
and  Measurement  of  Picosecond  Pul s  cwidth .  When  two  fundamental 
beams  have  properly  spatial  and  temporal  overlapping  inside  a  KDP 
crystal,  second  harmonic  signal  will  be  generated.  This  particular 
situation  is  well  described  previouslv  in  chapter  1^  To  ensure  the 
temporal  overlapping  inside  the  KDP  crystal^  Q-switched  Nd:glass 
was  used  as  a  source  of  excitation  since  it  had  relatively  long  pulse 
duration.  Q-switched  pulse  was  produced  by  means  oi  Pockels  cell 
(Korad  K-QS2  )  and  polarizer  stag  inside  the  laser  cavity.  The  arrange¬ 
ment  for  the  monitor  arm  was  exactly  the  same  as  described  in  the 
previous  setup  .  The  laser  was  again  regulated  by  a  slit  of  2  mm  wide 
and  5  mm  high.  A  schematic  of  the  experimental  arrangement  is  shown 
in  Fig.  3-5.  After  passing  through  the  slit  the  fundamental  beam  was 
then  reflected  by  a  mirror  having  anti  re  flection  coating  at  \  =  1.  06  pm. 

The  two  fundamental  beams  were  then  made  incident  upon  th  e  KDP  crystal.  The 
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KDP  crystal  and  the  detection  arm  for  SHI  from  the  KDP  were  fixed 
and  aligned  with  the  help  of  He-Ne  gas  laser  which  shone  tangentially 
through  the  surface  of  the  beam  splitter  (see  the  set  up  in  Fig.  3-5  ). 

The  KDP  crystal  is  oriented  in  such  a  way  that  its  optic  (z  )  axis  was 
parallel  to  the  entrance  surface  and  in  the  plane  of  reflection.  This 
orientation  of  the  KDP  would  yield  the  noncollinear  phase  matching 
condition,  which  was  experimentally  confirmed  and  will  be  discussed  in 
Chapter  IV.  Angle  of  incidence  9  of  each  fundamental  beam  was  equal 
and  regulated  by  the  rotation  of  reflecting  mirrors.  The  angular  turning 
platforms  on  which  both  reflecting  mirrors  were  situated  had  a 
144:1  reduction  gear  ratio  and  a  vernier  scale  permitting  a  variation  of 
91  in  step  of  1  minute  (1/60°).  The  polarizations  of  both  incident  beams 
were  in  the  direction  normal  to  the  plane  of  reflection.  In  another 
words  they  are  in  r  110]  direction  with  respect  to  KDP  crystallographic 
axes.  The  detection  arm  was  exactly  the  same  as  in  the  previous  setup, 
however,  it  was  fixed  in  a  normal  direction  with  respect  to  the  exit 
surface  of  the  KDP  crystal.  The  Q-switched  pulse  shape  of  the  fundamental 
beam  was  monitored  by  Tektronic  519  Oscilloscope  which  in  turn  provided 
triggering  signal  to  the  Tektronic  dual  beam  oscilloscope  on  which  both 
harmonic  signals  from  the  monitor  and  signal  channel  (from  KDP)  were 
displayed  and  photographed. 

The  experimental  arrangement  for  measuring  the  picosecond 
puls ewidth  was  almost  the  same  as  in  the  TBSM  experiment.  However, 

the  source  of  excitation  was  a  mode  locked  Nd:glass  laser  which  provided 
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picosecond  pulse  trains.  The  experimental  setup  is  shown  in  Fig.  3-6. 
Temporal  delay  for  each  beam  was  introduced  between  the  beam 
splitter  and  the  reflecting  mirror.  One  delay  was  fixed  and  the  other 
variable.  A  system  of  optical  delay  was  composed  of  three  right  angle 
prisms  as  shown  in  Fig.  2  -10.  Variable  temporal  delay  was  achieved  by 
mounting  the  prism  on  a  very  sensitive  translational  stage, 

which  could  translate  in  step  of  1  mil  (0.  0254  mm),  roughly  equal  to  1/10 
Picosecond.  The  KDP  crystal  with  its  orientation  used  in  the  investigation 
was  the  same  as  in  TFSM  case.  In  this  experiment  a  single  fundamental 
beam  could  not  produce  SHI  in  the  normal  direction  to  the  exit  face  of  the 
KDP  crystal.  It  required  both  spatial  and  temporal  overlapping  of  the 
two  beams  in  the  crystal.  Spatial  overlapping  was  easily  arranged  by  proper 
ly  rotating  the  two  reflecting  mirrors.  However  temporal  overlapping 
ould  be  °"ly  by  straight  forward  scanning  the  temporal  delay 

until  SHI  was  observed.  By  varying  temporal  delay  and  recording  the 

corresponding  SHI,  the  autocorrelation  function  was  mapped  out  and  the 
picosecond  pulsewidth  could  be  deduced. 
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CHAPTER  IV 

EXPERIMENTAL  RESULTS 

In  this  chapter,  experimental  results  of  second  harmonic  generation 
(SHG)  in  reflection  and  transmission  from  KDP  crystal  of  various 
crystallographic  orientations,  using  a  mode  locked  and  Q-switched 
Nd : glass  lasers,  are  presented.  Furthermore,  an  experimental 
technique  of  utilizing  the  polarization  properties  and  phase  matched 
condition  in  a  KDP  crystal  for  measurement  of  picosecond  pulsewidth 
is  reported  and  accounted  for.  The  results  of  the  present  work  are 
divided  into  four  major  parts  namely  the  SHG  in  reflection,  transmission, 
SHG  by  two  beam  spatial  mixing  (TBSM),  and  measurement  of  picosecond 
pulsewidth  of  Nd:glass  laser.  Those  results  are  presented  in  the 
following  paragraphs. 


A‘  Second  Harmonic  Generation  (SHG)  in  Reflection. 

Nonphase  Matchable  SHG  at  Total  Reflection.  The  KDP  crystal 
used  in  this  experiment  has  dimensions  25  x  15  x  8  mm.3  The  entrance 
surface  is  the  25  x  15  mm2  face.  The  face  normal  of  the  entrance 

surface  is  along  the  optic  axis  which  is  [  001]  direction.  The  total 

reflection  of  the  fundamental  laser  beam  and  for  harmonic  beam  are 
achieved  by  means  of  optically  denser  linear  fluid  1-bromonaphthalene  as 
previously  described.  The  polarization  of  the  fundamental  beam,  as 
indicated  in  Fig.  4-1,  is  along  the  f  110  ]  direction. 

The  reflected  second  harmonic  intensity  generated  from  the  crystal 
was  observed  as  a  function  of  the  angle  of  incidence  0.  for  the 
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crystallographic  orientation  shown  in  the  inset  of  Fig.  4-1.  The  KDP 
crystal  is  a  noncentrosymmetric  cubic  crystal  and  its  linear  susceptibility 
is  a  scalar.  The  theory  developed  by  Bloembergen  and  Pershan(12), 
which  is  given  in  Chapter  II,  can  be  applied  directly  here.  The  reflected 
harmonic  field  polarized  in  the  plane  of  reflection  is  given  by  equation 


(2.  33a) 


E  n  =  4tt  PNLS  F  NL 
R>n  R,ll 


(2.33a) 


where  the  nonlinear  Fresnel  factor  FRN^  is  given  by  equation  (2.35a) 


fNL  SinQS  SinQT  sin(*  +9s  +  V 

R,U  ®~R  Sin0R  sin^T+9R)c0s(9T‘9R)sin(9-+q-) 


T  S' 


(2.  35a) 


The  angles  9R>  9g  and  9^  are  related  to  9.  by  equation  (2-29) 
which  is  recast  into  the  forms 


Sin0R  =  r  %<“>/% 

sinqs  =  [n  (a,)/ncr(«,)]  sin9. 

sin9T  =  f nlicj(uu)/ncr (2au  )]  sin9. 

where  n^  stands  for  reflective  indicies  of  the  KDP  crystal  which ar 

in  chapter  III. 


(4.1) 


e  given 


Using  indicies  of  refraction  for  the  liquid  and  the  KDP,  The  critical 
angles  of  the  fundamental  and  second  harmonic  beams  can  be  found, 
according  to  equation  (4.1)  as 


(4.2) 


c  r 

For  9.>  9  (u» ),  9g  becomes  complex  and  one  should  use 
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cos  9g  =  if  (sin9. /sin9Cr(uj))2  -l]^2 


(4.  3) 


inequations  (2.32b),  (2.  35a)  and  (2.39).  Similarly,  for  9.  >  9Cr(2(JU )  , 

9^,  should  be  expressed  as 


cos9,r  =  ir(sin9./sin9Cr(2(JU))2  -l]1/2 


(4.4) 


(2.  39) 


According  to  the  theory,  the  reflected  second  harmonic  intensity 
!r(2uj)  ,  given  by  equation  (2.39),  is 

Ir(2»)=,c/8„,./Tr  |Eo|4dd'(4„X36NL,V(fRV  x 

lF  I  I  FpNL|2  cos9r(cos9.)"1 

The  solid  curve  drawn  in  Fig.  4-1  is  a  plot  of  equation  (2.  39)  on  a 

relative  scale  and  is  calculated  from  the  last  five  factors 

(f  L.2  |  L,  4  ,  NL ,  2  _i 

i  lF  I  I  fr  I  cos9R(cos9J 

This  is  appropriate  since  the  remaining  factors  in  equation  (2.39)  can 

be  treated  as  constants  in  the  experiment.  The  vertical  scale  in  Fig.  4-1 

is  adjusted  tothedata.  The  experimental  points  are  in  striking  agreement 

with  the  computed  theoretical  curve  which  displays  nonanalytical  singularities 

at  angle  of  incidence  9.  =  9Cr  (eu)  and  9.  =  9Cr(2w)  respectively.  These 

two  singularities  are  expected  since  at  these  angle,  cos9s  and  cos* 

respectively,  change  froma  real  value  to  a  pure  imaginary  value  .  We  would 

like  to  point  out  specifically  that  the  data  points  even  reproduce  two 

cusps  in  Fig.  4-1.  This  experimental  features  were  not  observed  in  the 

..  (28,29) 

earlier  measurements  of  this  sort  _  We  attribute  the 


ED 


70 

success  here  to  the  fact  that  in  our  experiments  the  mode  locked  laser 
beam  was  used.  The  unique  phase  relationship  among  the  axial  modes 
in  the  fundamental  beam  eliminated  most  of  the  intrinsic  quantum 
fluctuations  which  are  typical  of  nonlinear  optical  experiments  \ye 

thus  emphasize  once  more  the  importance  of  using  the  mode  locked  laser 
in  these  experiments. 

The  enhancement  of  the  reflected  intensity  arises  mainly  from  two 

sources.  First,  the  linear  Fresnel  factor  F  near  the  critical  angle  is 

larger  than  that  away  from  this  angle  by  about  a  factor  of  two.  This  will 

give  a  factor  of  sixteen  in  the  reflected  second  harmonic  intensity  I  (2(d). 

R 

Second,  the  nonlinear  Fresnel  factor  F  NL  defined  in  equation  (2.35a) 

K  t  11 

in  the  neighborhood  of  these  critical  points  as  dominated  by  the  term 

rsin(R  9g)]  1  and  it  is  larger  than  that  away  from  these  points  by 

(88 ) 

about  a  factor  of  three  .  Therefore,  after  all  factors  are  accounted  for, 

the  additional  enchancement  of  I  (2(d)  in  the  neighborhood  of  the  critical 

angles  will  be  about  two  orders  of  magnitude.  This  result  is  in  good 

agreement  to  the  theoretical  curve. 

The  physical  interpretation  of  the  enhancement  of  nonlinear  Fresnel 

factor  is  that  momentum  matching^  ’  ^  for  the  wave  propagating 

parallel  to  the  surface  inside  the  nonlinear  medium  is  important  in 

determining  the  reflected  second  harmonic  intensity  I  (2(d)  at  the 

critical  angle.  The  reflected  second  harmonic  intensity  near  the 

critical  angle  is  generated  by  the  polarization  in  a  surface  layer  with 

1/2 

thickness  of  about  (\f  )  where  f  .  is  defined  as  the  coherent 

coh  coh 
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The  enhancement  of  the  reflected  intensity  arises  mainly  from  two 

sources.  First,  the  linear  Fresnel  factor  FL  near  the  critical  angle  is 

larger  than  that  away  from  this  angle  by  about  a  factor  of  two.  This  will 

give  a  factor  of  sixteen  in  the  reflected  second  harmonic  intensity  I  (2uu). 

R 

Second,  the  nonlinear  Fresnel  factor  F  NL  defined  in  equation  (2.35a) 
m  the  neighborhood  of  these  critical  points  as  dominated  by  the  term 
s^n(^'p^  V]  and  it  is  larger  than  that  away  from  these  points  by 

ZOO) 

about  a  factor  of  three  .  Therefore,  after  all  factors  are  accounted  for, 
the  additional  enchancement  of  IR (2uj)  in  the  neighborhood  of  the  critical 
angles  will  be  about  two  orders  of  magnitude.  This  result  is  in  good 
agreement  to  the  theoretical  curve. 

The  physical  interpretation  of  the  enhancement  of  nonlinear  Fresnel 

factor  is  that  momentum  matching  ’  ^  for  the  wave  propagating 

paral  el  to  the  surface  inside  the  nonlinear  medium  is  important  in 

determining  the  reflected  second  harmonic  intensity  I  (2uu)  at  the 

R 

critical  angle.  The  reflected  second  harmonic  intensity  near  the 

critical  angle  is  generated  by  the  polarization  in  a  surface  layer  with 

1/2 

thickness  of  about  (\*coh)  where  l  ^  is  defined  as  the  coherent 
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ngth.  In  general,  the  reflected  second  harmonic  intensity  at  angles 

away  from  the  critical  angle  is  only  generated  from  a  layer  about 
/•/  e  thick. 

2'  — ase  Matched  Second  Harmonic  Generation  (SHG1  at  Total 
Rjn££ii2n.  The  KDP  crystal  used  in  this  case  has  dimensions  of 

25  *  'Z  *  18  mm  '  The  entrance  surface  is  25  x  12  mm2  face  and  it  has 
a  face  norma,  in  the  [111]  direction.  The  crys.a,  was  immersed  in 
•he  optically  denser  linear  fluid  1-bromonaphthalene.  The  special  cut 
was  specially  made  for  the  crystal  in  such  a  way  that  the  phase  matching 
direction  is  along  the  entrance  surface  of  the  crystal.  This  can  be  done 

by  cutting  the  crystal  such  that  the  optic  <z)  axis.  which  is  in  the  same 

direction  of  nonlinear  polari  zafinn  .pNLS  i 

polarization  P  ,  makes  a  phase  matching  angle 

’m  l°  th”  SUrfaCe'  The  an8le  was  computed  by  using  equation 
(2.  46)  and  indices  of  refraction  of  KDP  crystal.  It  was  found  that  at  the 

wavelength  of  interest  the  angle  ^  takes  value  of  41.  2°.  The  reflected 
second  harmonic  intensity  IR(2«,)  generated  by  the  crystal  was 

Observed  as  a  function  of  the  angle  of  incidence  ,.  for  the  crystallographi 

orientation  shown  in  the  inset  of  Fie  4-?  Up  r  i 

rig.  4-4.  The  fundamental  field  was 

polarized  along  the  flTo]  direction  and  propagated  as  an  ordinary  ray. 
According  to  equation  ,2.  31a)  the  nonvanishing  nonlinear  polarization 
was  along  the  optic  axis  and  therefore  the  harmonic  field  was  polarized 
in  the  plane  of  reflection  and  propagated  as  an  extraordina  ry  ray. 

The  eXperimen fa,  result  in  Pig.  4-2  shows  that  the  reflected  second 
h-rmonic  intensity  increases  by  about  three  orders  of  magnitude  when 
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the  angle  of  incidence  is  changed  by  a  few  tenth  of  a  degree  in  the  vicinity 
of  the  critical  angle.  There  is  one  maximum  (peak)  value  of  the  reflected 
harmonic  intensity  and  it  is  exactly  at  9*  =  0°r (m)  =  66.78°  in  agreement 
to  the  theoretical  prediction.  The  observed  maximum  reflected  intensity 
is  about  3  0  times  of  magnitude  larger  than  reflected  intensity  in  the 
nonphase  matching  case  described  in  the  previous  section.  This  is  due 
to  more  enhancement  of  the  reflected  harmonic  intensity  by  nonlinear 
Fresnel  factor  in  the  phase  matching  case. 

The  enhancement  can  be  mathematically  explained  in  the  view  of 
equation  (2.  39)  and  (2.  35a).  According  to  the  crystallographic  cut  of  the 
crystal,  the  value  of  the  nonlinear  Fresnel  factor  F^  ^  tends  to  infinity 
when  the  angle  of  incidence  9*  =  qCr(tlJ)  .  At  this  condition  we  have 

^  T  where  the  phase  matching  condition  prevails  and  thus 

sin(^s+9T)=0  in  equation  (2.  35a).  However,  in  practice,  there  are 
several  factors,  e.  g.  beam  divergence,  walk-off  effect,  that  limit  the 
reflected  harmonic  intensity  from  reaching  infinite  value.  Furthermore, 
the  singularity  causing  the  divergence  at  critical  total  reflection  predicted 
by  the  theory  for  this  situation  has  been  removed  by  another  treatment 
Riven  by  Shih  and  Bloembergen(91).  They  used  the  Green'  s  function 
technique  to  describe  the  generated  harmonic  field,  with  due  attention 
being  paid  to  the  effects  of  the  finite  beam  diameter.  According  to  the 
theory  based  on  this  technique,  the  second  harmonic  intensity  at  the 
critical  angle  QCr (w)  does  have  a  finite  value  and  it  is  in  agreement  to 
the  experimental  result  that  the  enhancement  of  the  reflected  second 


ha 


rmomc  intensity  is  three  to  four  orders  of  magnitude  when  the  angle 
of  incidence  is  changed  by  a  few  tenth  of  a  degree  in  the  vicinity  of 


c  r 

0  (t). 


The  theoretical  curve  in  Fig.  4-2  has  been  calculated  using  again 


equation  (2.  39)  with  =  q  and  ry  =  9  -  41  ?°  ,  a 

S  T  *  ’  4  ,Z  and  assuming 


ncr(,)  ne(2(V  )  1.4943  without  including  angular  variations  for  k 

-*  S 


and  kT.  This  is  justifiable  that  in  the  vicinity  of  the  critical  angle  the 


extraordinary  index  of  refraction  for  second  harmonic  frequency  is 
slowly  varying.  A  complete,  but  very  involved  expression,  has  been 

given  by  Fischer  (891  for  rigorous  treatment.  The  theoretical  curve  is 
again  calculated  from  the  last  five  factors. 


|fR  I  lF  I  I  FRj  „""l"  COs9r(cos9.) 


NL ,  2 


-1 


Of  equation  (2. 39).  The  vertical  scale  in  Fig.  4-2  was  adjusted  to  the 
data.  Note  the  striking  agreement  between  the  experimental  points  and 
eorefical  curve,  which  predicts  anomalously  high  reflected 
harmonic  intensity  ,M.  =  Furthermore,  i,  is  seen  that  the 

variation  in  the  immediate  neighborhood  of  the  critical  angle  is  rather 
well  described  by  equation  (2.  39)  for  the  geometry  used.  In  the  region 
where  ft1  is  greater  than  0Cr(„t) ,  and  also  cosd,.  and  cos9T  take 
pure  imaginary  values,  the  reflected  harmonic  intensity  is  still 


comfortably  detectable  and  it  tends  to  decrease  monotonously  to  aero 


value  when  Q  3 DDrosrli  Qf)^  prL: _  __  * 

i  Pproach  90  *  Thls  can  be  understood  physically  because 


bevond  the  critical  angle  there  is  still  evanescent  wave  of  the  fundamental 


am  wh.chw.ll  create  the  nonlinear  polarisation  PNLS(2«,)a„d will  radia te  back 
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into  the  linear  medium  the  reflected  harmonic  beam.  The  enhancement 
due  to  phase  matching  at  critical  angle  corresponding  to  this  geometry  is 
very  useful  for  studying  nonlinear  optical  properties  of  the  medium 
which  absorbs  at  the  fundamental  and/or  second  harmonic  frequencies. 
Since  all  the  information  will  be  analyized  via  reflected  Harmonic  intensity. 

To  facilitate  better  physical  understanding  in  connection  with  reflected 
second  harmonic  phenomena,  one  can  consider  the  drawn  curve  in 

4-2  as  a  limiting  case  of  that  in  Fig.  4-1.  The  separation  between 
the  two  cusps  in  Fig.  4-1  depends  upon  the  degree  of  phase  mismatching 
in  (he  neighborhood  of  the  critical  angle.  As  the  condition  for 
perfect  phase  matching  is  approached,  the  two  cusps  in  Fig.  4-1  will 
move  closer  towards  each  other  and  the  reflected  harmonic  intensity 
becomes  larger  until  finally  these  two  cusps  collapse  into  one  peak 
and  the  intensity  tends  to  the  highest  value  .  This  explanation 

was  also  carried  out  to  reflected  third  harmonic  generation  (THG)  by 

u  .  ,(30) 

Bey  et  al 

Another  experimental  investigation  in  relation  to  phase  matching  at 

the  critical  angle  was  performed.  Here  the  same  KDP  crystal  of  the 

same  crystallographic  cut  was  employed.  Howeverjthe  crystal  was 

rotated  180  about  the  face  normal  which  is  in  r  1 1 1  ]  direction  as  in  the 

previous  case.  The  reflected  harmonic  intensity  I  (2tu)  generated  from 

the  crystal  was  observed  as  a  function  of  the  angle  of  incidence  0. 

i 

The  crystallographic  orientation  and  the  polarization  of  the  fundamental 
beam  are  shown  in  the  inset  of  Fig.  4-3.  Again  in  this  case  the  fund- 
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amental  beam  propagated  as  ordinary  ray  and  the  reflected  harmonic 
beam  polarized  in  the  plane  of  reflection  was  extraordinary  ray. 

The  experimental  result  in  Fig.  4-3  shows  that  the  reflected 
harmonic  intensity  increases  by  about  two  orders  of  magnitude  when  the 
angle  of  incidence  changes  by  a  few  tenth  of  a  degree  from  the  critical 
angle.  The  maximum  reflected  harmonic  intensity  occurs  at  the  critical 
angle  as  predicted  by  the  theory.  The  theoretical  curve  in  Fig.  4-3 
was  calculated  in  the  same  way  as  in  Fig.  4-2.  The  vertical  scale  was 
adjusted  to  the  data.  Note  the  striking  agreement  between  the 


experimental  points  and  the  theoretical  curve,  which  predicts 
anomalously  high  reflected  harmonic  intensity  at  9.  =  9Cr(„)) 

When  the  result  from  Fig.  4-3  is  compared  to  that  from  Fig.  4-2, 
it  is  found  out  that  the  overall  shape  of  the  peak  in  Fig.  4-3  is  broader 
and  in  addition  Ihe  peak  intensity  s  lower  than  that  in  Fig.  4-2  by 
about  13  times.  This  discrepancy  can  be  well  understood  by  the  fact 


that  KDP  crystal  is  an  uniaxial  crystal  possessing  variable  extraordinary 


index  of  refraction. 


In  general,  one  may  expect  different  reflected  harmonic  intensities 
from  the  two  geometries  as  shown  in  Fig.  4-2  and  4-3,  respectively. 
These  two  geometries  are  identical  except  that  in  cne  case  the  KDP 
crystal  has  been  rotated  by  180°  about  its  face  normal  with  respect  to 
each  other.  The  two  situations  corresponding  to  Fig.  4-2  and  4-3  can  be 
shown  m  Fig.  4-4  where  the  normal  (index)  surfaces  of  ordinary  ray  and 
second  harmonic  extraordinary  ray  are  fully  indicated  for  both  cases. 


This  is  equivalent  to  an  illumination  by  a  beam  of  incident  from  a 
symmetric  position  on  the  other  side  of  the  face  normal.  In  Fig.  4-4, 
where  the  drawing  indicates  the  conservation  of  the  tangential  components  of 
lc  vectors.  The  k  vectors  with  subscripts  a  and  b  correspond  to 
the  geometries  in  the  insets  of  Fig.  4-2  and  4-3,  respectively.  At 

a  -*  , 

total  reflection,  kq  =  OA  ,  kg  =OC.  According  to  the  normal 
(index)  surfaces  and  boundary  conditions  of  k>  s  vectors,  it  is  required 
that  kT  may  take  two  values,  i.e.  ]<T  =  k  3  =  OA  or  k  3  =  OB  . 

However,  the  former  solution  of  k^  =  0A  would  give  a  ray  velocity 
propagating  out  of  the  nonlinear  medium;  i.  e.  ,  the  energy  propagates 
out  of  the  medium.  This  is  not  a  physically  allowed  solution.  Thus  the 

only  physically  allowed  solution  is  t k^  =  O  B  and  k  b  =  O  C.  They 
correspond  to  the  incident  k,3  and  k.b  waves,  respectively. 

In  the  formal  treatment  of  the  theory,  Armstrong  et  al(11)  and 
Klemman  have  represented  the  laser  and  second  harmonic  light  by 
unbounded  plane  waves.  Armstrong  et  al(H)  pointed  out  that  the  inter¬ 
action  between  light  waves  of  finite  aperture  takes  place  along  the 
direction  of  energy  flow.  In  his  paper,  Kleinman 04 ]  also  showed  that 
PHG  will  ultimately  be  limited  by  slightly  different  directions  of 
propagation  of  the  energy  of  the  laser  and  second  harmonic  beams  .  The 
effect  is  called  double  refraction  in  a  uniaxial  crystals.  In  fact,  the 
light  wave  in  the  nonlinear  medium  KDP  in  the  two  geometries  have  very 
narrow  apertures  near  critical  angle.  Therefore,  to  examine  the 
phenomena  of  SHG  at  phase  matched  condition  in  a  uniaxial  crystal,  e.  g. 


KDP,  it  is  important  to  know  the  degree  of  overlapping  between  the 
di  rection  .of  energy  propagation  of  polarization  wave  k  and  radiation 
wave  kT>  It  has  been  shown  by  Boyd  et  al(35)  that  for  any  direction  of 
phase  propagation  the  direction  of  energy  propagation  (Poynting  vector) 
is  given  by  the  normal  vector  to  the  normal  (index)  surface.  For  the 
ordinary  wave  the  propagation  vector  and  the  Poyting  vector  are  parallel. 
However,  for  the  extraordinary  wave  the  Poynting  vector  deviates  fro, 
the  phase  propagation  direction  by  an  angle  p  ,  called  the  angle  of 
double  refraction,  the  angle  p  is  given  by^35^ 

1  e  .  r  1  1 


>m 


1  e  r  1 

tan  o  =  7  n  (9)  i - 

£  c<\)  t  _  e 


r"2>/2>l2  r"°J2  ' 


1„. 


sin29 


(4.  5) 


where  ‘  is  the  angle  between  the  phase  propagation  direction  and  the 

optic  (z)  axis.  Furthermore,  the  angle  Cp  in  Fig.  4-4  can  be  obtained 
from 
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rn 6  (P)l2  rn_°  ]2 

2uu  2<ju  J 


cos9  + 


•  2n 

sin  9 


(4 .  6a ) 


n  (41.  2°-co)  coscp  =  n° 

^  'JU  (j) 

From  (4.  6a,  b)  the  angle  cp  and  as  a  cons  equence  nG  (41.2°-cp) 

2  m 


(4.  6b) 


were 


found  as 


cp  =  1.12 


n2uj  (4L2°-1-  12°)  =  nG  (40.  08°)=  1.49465 


2uu 


The  angle  0,  was  found,  upon  substitution  nG  (40  08°)  and 

2(u 

appropriated  value  of  indices  of  refraction  of  KDP,  as 


(4.  7) 


R1 


Oj  =  1-  07° 

In  the  same  way  was  calculated  as 
o2  =  1.08° 


(4.  8) 


(4.  9) 


Therefore, the  angle  between  the  directions  of  energy  propagation  of 
and  1<T  corresponding  to  the  geometry  of  Fig.  4-2  is 

CP  -  Oj  =  1.  12°  -  1.  07°  =  0.  05° 

whereas  the  angle  between  the  direction  of  energy  propagation  of 

-*b  S 

and  kT  corresponding  to  the  geometry  of  Fig.  4-3,  is  o2  =  1.  08°. 

Thus  in  the  former  case,  the  polarization  and  radiation  waves  will  have  a 

larger  interaction  volume  and  hence  will  produce  more  intense  second 
harmonic  wave. 

This  qualitative  analysis  via  double  refraction  phenomena  explains 
the  experimental  results  shown  in  Fig.  4-2  and  4-3.  However, the 
double  refraction  effect  corresponding  to  the  idealized  geometry  of  phase 
matching  shown  in  Fig.  4-4  probably  cannot  account  for  the  experimental 
factor  of  13  times  increase.  The  effect  mentioned  above  depends  very 
much  on  the  actual  value  of  the  angle  cp  and  the  sizes  of  the  fundamental 
and  harmonic  beams  near  the  critical  angle  inside  the  nonlinear  medium. 

Another  possible  explanation  of  the  different  intensities  when  the 
KDP  crystal  is  turned  by  180°  can  be  based  upon  the  beating  of  different 
ordinary  ray  directions  in  the  fundamental  wave.  This  idea  was  first 
given  by  Maker  et  al(34).  In  this  explanation,  we  assume  that  the  crystal 
was  not  cut  precisely  that  the  phase  matching  direction  is  not  along  its 
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surface  and  in  addition  the  fundamental  beam  is  allowed  to  have  a  slight 
divergence.  The  illustration  is  shown  in  Fig.  4-5  where  symbols  of 
subscripts  a  and  b  denote  the  same  geometries  as  previously 
described  in  the  early  analysis.  According  to  this  case,  the  experimental 
maximum  intensity  is  not  obtained  for  the  idealized  geometry  of  Fig.  4-4, 
however,  the  maximum  occurs  when  the  phase  matching  direction  0  A 


make  a  small  angle  5  with  the  surface  as  shown  in  Fig.  4-5.  The 


nonlinear  polarization  is  created  from  the  mixing  of  slightly  divergent 


~*a  '  -*a" 

pairs  ol  fundamental  rays  with  wave  vectors  k  and  k  .  The  allowed 

S  S 


perfect  phase  matching  condition  is 


’  /  -*  // 
a  .  a 


:T  =  kS  +  kS 


This  condition  can  be  satisfied  if  5,  >9  ,  as  discussed  bv  Maker 

k  «■>  1 


,(34) 


T 


m 


et  al  ,  where  5^  is  the  angle  between  k^  and  the  optic  (z  )  axis 


T 


When  the  crystal  is  turned  by  180  around  its  face  normal,  the  wave 


V 


vectors  of  the  incident  fundamental  rays  are  represented  by  k  and 

%"  S 

k<-  •  For  this  situation  the  phase  matched  condition  is  not  satisfied  and 


as  a  consequence  the  reflected  harmonic  intensity  will  be  lower.  The 
analysis  based  on  the  discussion  by  Maker  et  al^4'  seems  to  agree  with 
the  experimental  data  obtained  in  the  case  of  nonphase  matched  and 
phase  matched  conditions  described  in  Fig.  4-1  and  4-3,  respectively. 

It  was  experimentally  found  that  ratio  of  the  harmonic  intensities  at  the  critical 
angle  for  phase  matched  to  nonphase  matched  conditions  is  about  30  times, 
which  is  in  the  same  order  of  magnitude  (13  times)  as  obtained  upon 
comoaring  Fig.  4-3  to  Fig.  4-2 


\ 
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FIGURE  4  -  5 


3.  Nonlinear  Brewster1  s  Angle.  The  agreement  of  experimental 

data  with  theory  for  those  separate  experiments  described  in  the  previous 

sections  encourages  verification  of  a  more  subtle  nonlinear  optical 

postulate.  In  particular,  we  are  led  to  ask  if  a  nonlinear  Brewster'  s  angle 

for  the  nonlinear  transparent  medium  can  be  found  experimentally. 

Even  though  the  nonlinear  Brewster'  s  angle  for  absorbing  medium,  e.g.  , 

(27) 

GaAs,was  first  confirmed  by  Chang  and  Bloembergen  ,  it  is  still 

important  to  carry  the  postulate  to  the  transparent  regime.  An  uniaxial 

KDP  crystal  is  known  to  be  transparent  for  both  fundamental  and 

second  harmonic  frequencies  of  Nd:glass  laser.  Besides  its  transparent 

property,  it  is  also  an  ideal  nonlinear  material  in  the  sense  that  it 

crystallographic,  point  group  properties  and  nonlinear  optical  characteristics 

concerning  SHG  are  well  understood.  The  KDP  crystal  used  in  the 

investigation  is  the  same  as  one  used  in  the  phase  matching  condition  at 

total  reflection  described  in  Fig.  4-3.  Again  the  crystal  is  immersed  in 

the  optically  denser  linear  fluid  1-bromonaphthalene.  The  use  of  the 

liquid  is  to  facilitate  the  achievement  of  observation  of  Nonlinear 

Brewster'  s  angle  since  without  optically  denser  linear  medium  from 

which  the  fundamental  is  incident,  the  condition  for  Nonlinear  Berwster'  s 

angle  will  never  be  achieved.  The  face  normal  of  the  entrance  surface 

of  the  crystal  is  in  Till]  direction  and  the  fundamental  wave  is  polarized 

along  rl  1  0  ]  dir ection.  According  to  (2. 31a)  the  nonlinear  polarization 

P  will  be  in  r001]  direction  or  optic  (z )  axis.  The  reflected 

harmonic  intensity  I  (2 oj )  generated  from  the  crystal  was  observed  as 

R 
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a  function  of  the  angle  of  incidence,  which  varies  from  20°  to  75°.  The 
experimental  result  is  shown  in  Fig.  4-6.  In  this  particular  geometry 
of  the  KDP  crystal,  we  have  the  nonlinear  polarization  in  the  plane  of 
reflection.  The  electric  field  of  the  reflected  harmonic  wave  is  given 
by  equation  (2.  25)  .  This  equation  reveals  the  existence  of  a  Nonlinear 

p 

Bi  ewster1  s  angle  when  (2ui)  =  0  .  Or  in  another  words,  we  have 

.  _NLS  .  „  .2 

R  4nP  sin)  sin  9  smM  9  +  9  ) 

E  (2'jj)  =  0  =  - — - ; - — - 1 - 1 - § _ 

11  eRSlnqRsln(yT+  9s)sin(9T+6R)cos(9T-9R)  1U 

The  equation  (4.10)  is  true,  if 


sin(o-  +  9  +  0  )  =  0 

or  O'  +  9s+9T=0,n  (4.  11) 

From  the  Fig.  4-6,  one  has  a  +  9^  ,  which  is  the  angle  between  the 

nonlinear  polarization  and  the  face  normal  direction  inside  the  crystal, 

equal  to  -  48.  78  .  Therefo re  ,  from  equation  (4.  11 )  using  the  first 

condition  which  provides  the  realizable  value  of  9  we  have 

T 

9  =  48.  78° 

T 

Having  known  the  value  of  9  ^  =  48.  78  ,  the  corresponding  index  of 
refraction  of  an  extraordinary  ray  n®  (9  =  48.  78°)  is  found  via 

equation  (4.  6a)  to  be 

n2uu  (9T  =  48<  7SO)  =  L  47125 

Thus  the  Nonlinear  Brewster’  s  angle  can  be  found  via  nonlinear 
Snell  law  (2.  29)  as 


OF  INCIDENCE 


This  angle,  of  course,  is  nonphysical. 

Note  the  striking  agreement  between  experimental  result  and  the 

theoretical  prediction.  One  can  see  that  in  the  vicinity  of  01  =  42. 78°, 

there  exists  a  pronounced  dip  of  the  reflected  harmonic  intensity  I  (2<u). 

R 

This  evidently  confirms  the  existence  of  Nonlinear  Brewster'  s  angle. 

This  is  for  the  first  time  that  Nonlinear  Brewster1  s  angle  for  transparent 
medium  has  been  experimentally  demonstrated.  Notice  that  the  reflected 
harmonic  signal  in  the  neighborhood  of  Nonlinear  Brewster  angle  is  very 
low  level  and  at  the  angle  it  is  expected  no  IR  (2uu )  generated  out  of  the 
crystal.  By  increasing  the  incident  power  density  via  a  biconvex  lense  and 
using  the  low  level  signal  detection  method  in  which  the  reflected  harmonic 
signal  was  collected  and  averaged  over  25  to  50  laser  firings,  the 
pronounced  dip  at  9.  =  42.  83°  was  achieved.  Other  experimental  points 
agree  well  with  the  theoretical  curve  which  was  calculated,  as  in  previous 
case,  from  the  last  five  factors  of  equation  (2.39).  In  this  case  the 
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(he  angular  variation  of  n  (9  )  was  taken  into  consideration  in  computing 

the  theoretical  curve  since  we  are  now  dealing  with  much  larger  angular 

range.  From  Fig.  4-6,  again  the  anomalously  high  reflected  harmonic 

intensity  I  (2mi)  at  critical  angle  and  overall  agreement  between  data 

points  and  theory  from  2  0°<  q1  <  75°  confirm  the  nonlinear  optical  theory 

(12 ) 

developed  by  Bloembergen  and  Pershan 

The  physical  interpretation  of  the  existence  of  a  Nonlinear  Brewster'  s 
angle  can  be  understood  in  terms  of  clas  sical  dipole  radiation.  When 
the  fundamental  beam  is  incident  on  the  KDP  crystal  at  the  Nonlinear 

“♦‘NTT  q 

Brewster'  s  angle,  it  will  setup  a  nonlinear  polarization  P  in 
dii  e^tion  parallel  to  the  direction  of  propagation  of  the  reflected  harmonic 
wave  inside  the  nonlinear  medium.  According  to  the  classical  dipole 
radiation  theory,  there  would  be  no  radiation  seen  in  this  direction.  This 
nonradiating  wave  upon  refraction  back  into  the  linear  medium  would 
otherwise  give  rise  to  the  reflected  ray  in  the  direction  of  9  . 

The  experimental  points  and  the  corresponding  theoretical  curve  in 
the  vicinity  of  Nonlinear  Brewster'  s  angle  havea  pronounced  dip, in  contrast 
to  the  case  of  GaA»s  irradiated  by  ruby  laser,  performed  by  Chang  et  al^ 

This  can  be  explained  that  for  the  KDP  case,  the  crystal  has  real  and  small 
value  of  S (fJ) )  and  s(2(ji)  in  contrast  to  the  case  of  GaAs  which  has  large 
and  complex  values  ofe(ou)  and  e(2uu)  at  the  ruby  laser  line.  Maxwell'  s  equa 
tions  and  the  boundary  conditions  derived  for  a  transparent  medium  can  all 
be  extended  to  an  absorbing  medium^provided  that  the  dielectric  constants 
are  replaced  by  their  complex  values.  Then,  in  that  case, one  can  see 
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that  the  condition  for  (2(jj)  =  0  can  never  be  satisfied,  as 

1  /  z 

sinT  =  re  (2'u ) / G (|,! ) ^  sin  48.  78  can  never  be  achieved  physically 

if  s(2f|))  is  complex.  This  fact  was  confirmed  when  a  Ndrglass  laser  was 

used  instead  of  a  ruby  laser,  and  the  pronounced  dip  at'  the  Nonlinear  Brewster' 

(34) 

angle  for  GaAs  was  observed  by  Chang  et  al  .  This  is  because  the  second 
harmonic  emission  of  the  Ndrglass  laser  is  absoi bed  much  less  by 

GaAs  since  'he  imaginary  part  of  its  linear  dielectric  constant  is  1.  50. 

The  conclusion  to  be  drawn,  in  relation  to  the  reflected  intensity  of  the 
harmonic  wave  from  a  nonlinear  medium  KDP  crvstal  immersed  in  a 
op'ically  denser  linear  fluid  l-bromonaphthalene  ,  , s  that  the  observed 

angular  dependence  of  the  reflected  harmonic  intensity  and  ,  in 
particular, the  existence  of  Nonlinear  Brewster1  s  angle  of  a  transparent 
medium  are  in  good  agreement  with  the  theory  of  Bloembergen  and 
Pershan^  \  The  regime  of  an  evanescent  fundamental  wave  and 
reflected  second  harmonic  wave  have  been  demonstrated. 


B.  Second  Harmonic  Generation  (SHG)  in  Transmission. 

1.  Homogeneous  and  Inhomogeneous  Second  Harmonic  Waves  in  the 

Neighborhood  of  Critical  Angles.  Having  obtained  a  good  agreement 

between  the  experimental  data  and  theoretical  prediction  in  the  case  of  the 

nonphase  matchable  reflected  SHI  described  in  section  A-l,  it  was 

encouraging  to  investigate  the  SHG  in  transmission  near  the  critical 
c  r  c  r 

angles  5  (Zni)  and  9  (uj)  ,  respectively.  The  result  obtained  in 

this  investigation,  according  to  the  theoretical  prediction,  will 
complement  that  of  section  A-l.  The  gemetrical  situation  just 
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before  total  reflection  occurs  was  shown  in  Fig.  2-6  in  chapter  II.  The 
primary  fundamental  beam  is  transmitted  almost  parallel  to  the  surface  in 

/jo  \ 

the  nonlinear  KDP  crystal.  There  are  two  transmitted  harmonic  beams 
The  driven  polarization  wave  propagates  in  the  same  direction  as  the 

*4  -4 

transmitted  laser  beam.  It  has  a  wave  vector  k  =  2k  (ou)  and 

S  Laser 

represents  the  particular  solution  of  the  inhomogeneous  wave  equation. 

In  addition, there  i s  a  homogeneous  solution  with  wave  vector  k ^ (2 (ij )  . 

In  the  present  study  the  KDP  crystal  has  a  right  angular  corner  and  the  two 

transmitted  harmonic  beams  are  spatially  distinct  and  readily  observed 

separately.  According  to  the  particular  crystallographic  cut,  which  will 

be  described  later,  it  is  clear  that  the  beam  with  wa^e  vector  k  will 

T 

c  r  -4 

disappear  first  at  9^  =  9  (2oj )  und  the  ray  with  wave  vector  k^  will 

disappear  at  the  same  time  as  the  transmitted  fundamental  since  9Cr(oj) 

is  found  to  be  greater  than  9  (2uu)  .  A  s  the  angle  of  incidence  9  becomes 

i 

c  r 

larger  than  9  (ju  )  the  inhomogeneous  wave  will  disappear. 

Q 

he  KDP  crystal  used  has  dimensions  of  25  x  15  x  8  mm  .  Its  face 
normal  is  in  the  T001]  direction  (optic  axis).  The  homogeneous  and 
inhomogeneous  transmitted  harmonic  intens'-ty  were  observed  separately 
as  a  function  of  angle  of  incidence.  The  polarization  of  the  fundamental 
beam  and  the  orientation  of  the  crystal  are  shown  in  the  inset  of  Fig.  4-7. 

The  homogeneous  and  inhomogeneous  harmonic  intensities  are  given 
by  equation  (2.39),  respectively,  as 

V2<»>  =  (c/8ttK/Tr  lEo|4dd'(4TTX3^)2r,2(fb  |FL  |4  x 


Ft|  cos9t(cos9.) 


-1 


(4.12) 


TRANSMITTED  SHI  (RELATIVE  UNITS 


FIGURE  4-7 
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Is(2m)=  <c/8„K/eR|Eo|4  dd'  (4TTX36NL)2r)2(f^)|FL|4 

,  .2  "I 

|Fgl  cos  9g(co  s9. )  (4.13) 

The  drawn  solid  curves  in  Fig.  4-7  are  plots  of  equations  (4.  12)  and 
(4.  13)  for  homogeneous  and  inhomogeneous  SHI  and  are  calculated  from 
the  last  five  factors,  respectively, 

, ,  L  .2  |  L.4  ,  NL  ,2  -1 

h, s  |F  I  Ift.sI  coseT, 

The  linear  vertical  scale  in  Fig.  4-7  is  adjusted  to  the  data  as 
appropriat  for  each  case.  The  experimental  points  are  in  good  agree¬ 
ment  with  the  computed  theoretical  curve.  They  arr  in  excellent 
agreement  to  the  theoretical  prediction  that  the  homogeneous  and 
inhomogeneous  SHI  will  be  terminated  at  9°r (2'«)  and  0Cr((u)  , 
respectively.  The  critical  angles  ^Cr(2ijj)  and  0Cr(iu)  from  the 
experiment  are  the  same  as  found  in  the  previous  cases.  The  behavior 
of  the  two  SHI  near  the  critical  angles  as  shown  in  Fig.  4-7  can  be 
explained  in  the  following.  The  two  transmitted  SHI  are  the  net  result 
of  the  competing  effects  of  the  increasing  linear  Fresnel  factor  to  the 
fourth  power,  the  decreasing  nonlinear  Fresnel  factor  squared,  and  the 
rapidly  decreasing  cross  section  of  the  fundamental  inside  the  medium 
near  the  total  reflection.  As  to  these  competing  factors  the  homogeneous 

and  inhomogeneous  harmonic  intensities  fall  down  rapidly  toward 
c  r  cr 

9  (2(ii)  and  9  (tu)  .  The  reason  for  the  homogeneous  and  inhomogeneous 

c  r  cr 

SHI  terminated  at  0  (2(d)  and  9  (id)  respectively  is  that  when  9.  is 

greater  than  0C  (2uj)  and  9 ° r (UJ )  the  value  of  cosQ  and  cos9  become 

T  S 
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pure  imaginary  respectively.  As  a  consequence  the  two  transmitted 

harmonic  intensities  will  become  imaginary  which  are  not  physically 

allowable.  Furthermore,  one  can  see  from  Fig.  4-7  that  the  homogeneous 

SHI  disappears  at  9  (2uj)  while  inhomogeneous  SHI  still  persists  for 

2.  02  more  and  finally  is  terminated  at  0Cr(io)  .  In  the  region  between 

9  (2m)  s:  q.  ^  9  ((«)  which  is  2.  02°  interval,  one  can  comfortably 

observe  inhomogeneous  SHI.  One  can  apply  the  result  from  this 

investigation,  particularly  in  the  region  8 ° r (2c« )  5  0  <;  8 C r (^j)  to  other 

nonlinear  crystal  which  is  appropriately  cut.  Since  in  this  region  one 

will  obtain  only  inhomogeneous  SHI  which  has  direct  association  with  the 

nonlinear  polarization  source  term  P  and  the  nonlinear  susceptibility 
NL 

X  •  The  knowledge  of  inhomogeneous  SHI  in  this  particular  region  will 
directly  facilitate  the  study  of  PNLS  and  of  a  nonlinear  medium. 

2-  phase  Matched  Second  Harmonic  Generation  (SHG)  in  Transmission. 

SHG  under  phase  matched  conditions  and  in  a  normally  incident 
direction  was  performed  by  several  authors*33,  34\  However  in  the 
present  work  the  phase  matched  direction  is  no  longer  in  the  normal  incident 
direction.  The  SHG  in  transmission  under  this  phase  matched  condition, and 
away  from  it, were  experimentally  investigated  and  compared  to  the  prediction 
of  the  theory.  This  will  verity  the  theoretical  prediction  and  check  the 
theory  for  other  conditions  aside  from  previous  confirmations. 

The  KDP  crystal  used  in  the  experiment  was  the  same  one  as  in 
section  B-l.  The  tctaJ  transmitted  second  harmonic  intensity  due  to  the 
sura  of  homogeneous  and  inhomogeneous  transmitted  harmonic  intensities 
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is  observed  as  a  function  of  angle  of  incidence  which  varied  from  -15°  to 

50  .  The  fundamental  beam  is  polarized  in  the  f  110]  direction  and  the 

KDP  crystallographic  orientation  are  shown  in  Fig.  4-8.  The  phase 

matching  angle  9^  according  to  this  geometry  can  be  computed  by  using 

equation  (2.46)  and  indices  of  refraction  of  KDP  crystal  given  in  Chapter  III. 

The  phas  e  mat  ched  direction  was  found  to  be  9  =  41.  2  away  from  the 

crystal  face  normal  direction  (optic  axis).  Thus  the  angle  of  incidence 

cor  res  ponding  to  this  direction  is  given  by 

,  n®  (9  ) 

n  _  ■  f  2ii)  m  .  „  1 

9.  -  sin  i  - — - -  •  sin9  ' 

1  *  n  (<u)  mJ 

iiq 

where  n®  (9  )  =  n°  =  1.4943 

Zuu  m  uj 

r,  A  .  -If  1.4943  .  o ") 

Kence  9.  -  sin  •  sm41.2  [ 

9.  =  37.  27° 
i 

The  theoretical  curve  for  total  transmitted  second  harmonic  Intensity 
is  the  sum  of  homogeneous  and  inhomogeneous  intensities  given  by  (4.12) 
and  (4.13),  respectively 

^otal^  =  IS(2t',)  *  V2UJ)  <4-14) 

The  thereotical  curve  shown  in  Fig.  4-8  was  computed  from  the  sum  of 
the  last  five  factors  in  equation  (4.12)  and  (4.13)  and  it  can  be  mathematically 


expressed  as 


ITotal(2tW)^  (fT  }  IfL|4(cosV  ^  ,FT,n,2c°SV  lFS,Ill|2coses}  (4'15) 
The  expression  for  ITotal  =  +  IT  given  by  equation  (4.  14)  is 

justified  for  the  present  investigation.  For  this  geometry,  the  transmitted 
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intennties  were  measured  by  observing  the  transmitted  harmonic  waves 

emerging  from  the  back  face  of  the  crystal.  Since  the  two  harmonic 

beams  now  overlap  and  interfere, the  harmonic  intensity  is  a  function  of 

the  path  length  of  the  beam  inside  the  crystal.  If  the  entrance  and  exit 

faces  are  not  parallel  over  the  entire  cross  sectional  area  of  the  beam, 

the  interference  must  be  integrated  over  this  area.  The  result  of  this 

integration  gives  three  terms.  The  first  two  terms  are  exactly  the 

separate  homogeneous  and  inhomogeneous  intensities  given  by  (4.12)  and 

(4.  13)  respectively.  The  third  term  retains  an  oscillatory  dependence  of 

the  length  of  crystal  and  has  an  amplitude  determined  by  the  angle  of 

parallelism  of  the  two  faces.  If  the  two  faces  are  exactly  parallel, the 

magnitude  of  the  third  term  equals  to  the  sum  of  the  first  two  terms^°l 

If  the  path  lengths  of  the  two  harmonic  beam  differed  by  several  coherence 

lengths,  the  amplitude  of  the  third  term  is  negligible  compare  to  the  first 

two.  Since  the  KDP  crystal  used  in  the  experiment  is  an  uniaxial  crystal 

in  which  a  double  refraction  phenomena  is  intrinsic.  Therefore  the 

magnitude  of  the  third  term  corresponding  to  the  interference  of  the  two 

harmonic  beams  is  negligible  and  thus  the  expression  of  I  (2uj) 

To  tal 

given  by  (4.  14)  is  valid 

According  to  the  theory  for  general  crystallographic  orientation,  it 
is  anticipated  that  the  total  harmonic  intensity  will  exhibit  the  oscillatory 
feature  in  the  heighborhood  of  very  small  angle  of  incidence.  However, 
this  effect  has  not  been  observed.  Instead, the  total  harmonic  intensity 
near  normal  incidence  is  falling  rather  monotonically  toward  zero  value 
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at  Q  =  0  .  This  behavior  can  be  understood  in  terms  of  two  physical 
reasons.  First,  for  the  region  of  small  angle  of  incidence,  the  double 
refraction  still  persists  and  it  reduces  the  effectiveness  of  the  interference 
between  the  two  harmonic  beams  as  described  in  the  previous  paragraph. 
Second,  according  to  particular  crystallographic  orientation  of  the 
crystal  used  in  this  experiment,  it  is  noted  that  the  pNLS  {s  in  the 
same  direccion  of  the  crystal  face  normal.  In  terms  of  the  dipole  radiation 
point  of  view,  the  nonlinear  polarization  P  cannot  radiate  harmonic 
waves  in  the  direction  of  its  oscillation.  Therefore  it  is  anticipated  that 
no  harmonics  intensities  are  observed  in  the  direction  of  the  crystal 
face  normal  i.  e. ,  9.  =  0  .  Furthermore,  the  total  harmonic  intensity 

will  be  small  in  the  neighborhood  of  9.  =  0°  .  Thus  the  trend  of 
monotonic  decrease  of  the  total  harmonic  intensity  dominates  the 
anticipated  oscillatory  behavior  in  this  region. 

The  total  harmonic  intensity,  collected  by  means  of  the  technique 
described  in  details  in  chapter  III,  is  observed  as  a  function  of  the  angle 
of  incidence  which  varies  from  -15  to  55  .  Note  the  striking  agreement 
between  experimental  data  points  and  computed  theoretical  curve  of 
ITotal^tU  )  given  by  (4.14).  In  particular,  the  theoretical  prediction  con¬ 
cerning  the  phase  matching  angle  9.  =  37.27  and  zero  total  harmonic 
intensity  at  9.  =  0  are  confirmed.  In  addition,  oy  using  a  KDP  crystal 
and  a  mode  locked  Nd:glass  laser,  the  dynamical  range  of  the  dectable 
harmonic  intensity  is  of  about  10  orders  of  magnitude. 

Another  investigation  of  phase  matched  SHG  in  transmission  was 
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performed  by  using  a  KDP  crystal  of  different  crystallographic  cut.  The 

KDP  used  in  this  investigationha s  an  optic  axis  ,  which  is  r  00 1  ]  directed 

making  an  angle  41.2°  to  the  crystal  entrance  surface.  The  normal  to  the 

crystal  face  is  in  r  111]  direction.  The  crystal  is  immersed  as  usual  in  a 

linear  optically  denser  liquid  1-bromonapththalene.  The  fundamental  beam  is 

polarized  in  the  r  110  ]  direction.  The  total  harmonic  intensity  I  ^  ^ ( 2 ut ) 

was  observed  as  a  function  of  angle  of  incidence  which  varied  from  0°  to 

52.  5°.  The  crystallographic  orientation  of  the  crystal  and  the  result 

obtained  from  the  investigation  is  depicted  in  Fig.  4-9.  ,  The  theoretical 

curve  of  total  harmonic  intensity  I  ,  (2 uj  )  is  again  given  by  equation 

Total 

(4-14)  on  the  same  assumption  explained  in  the  previous  section.  The 

theoretical  curve  drawn  in  Fig.  4-9  was  computed  from  equation  (4-15) 

and  the  vertical  scale  is  adjusted  to  the  experimental  data. 

The  angle  of  phase  matching  ^  could  be  found  in  the  usual  way  by 

using  equation  (2.46)  and  indices  of  refraction  of  KDP  given  in  Chapter  III. 

The  phase  matching  angle  was  found  to  be  d  =  41.2°.  According  to  the 

m 

present  geometry,  there  exist  two  possible  phase  matched  directions.  The 
first  direction  is  along  the  entrance  surface  of  the  crystal.  This  direction 
is  associated  with  phase  matched  SHG  at  total  reflection  and  has  been 
investigated  in  section  A -2  previously.  The  second  direction  which  is 
still  41.  2°  away  from  the  other  side  of  optic  (z )  axis.  It  makes  an  angle 
9  ,  inside  the  KDP  crystal,  to  the  direction  of  the  face  normal.  The 
angle  9  is  given  by 

9  =  90°  -  (41.  2°+  41.  2°  )  =  7.  6° 
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The  corresponding  angle  of  incidence  9.  is  given  by 


.-ir'W  .  ! 

9  '  sm  IFTm — •  siii9t1 


liq 


where  n  (9  )  =  n  =  1.  4943 

4uu  m  a) 


Thus 


_  _  .  -1  f  1-  4943  ( 

i  "  S]n  076260  •  sin(7-6 


9.  =  6. 95 
i 


The  rapid  increase  of  the  total  harmonic  intensity  at  the  phase 


matching  angle  0^  is  due  to  an  enhancement  of  the  nonlinear  Fresnel 


,  *  „NL  NL 

aC°1S  S,  11  and  Ft,11  ?iven  by  (2.  35b)  and  (2.  35c)  for  polarization  and 


radiative  harmonic  waves  respectively.  This  rapidly  monoton 

increase  of  the  total  harmonic  intensity  dominates  the  oscillatory 

pattern  in  the  region  of  small  angle  of  incidence  so  that  no  oscillatory 

pattern  of  interference  between  the  homogeneous  and  inhomogeneous 

harmonic  beam  has  been  observed.  The  experimental  data  points 

are  in  striking  agreement  to  the  computed  theoretical  curve  shown  in 

Fig.  4-9.  It  is  noted  that  the  prominent  dip  of  the  curve  occurs  in  the 

neighborhood  of  9.  =  44°  as  expected  in  theory.  According  to  classical 

dipole  radiation,  the  inhomogeneous  or  homogeneous  harmonic  intensity 

will  become  minimum  when  the  propagating  vector  k  or  k  parallel 

o  T 


to  the  direction  of  PNLS  which  is  in  the  r001]  direction.  This  situation 


occurs  when  9  or  9  is  (90°  -  41  fl°  on. 

S  oT  ib  (vu  41.  «£  1-48.8  .  The  corresponding  angle 


of  incidence  9  is  thus  given  by 


MM - 


mMi  - iMi  turn  mui 


.  .  ■■ - - - - 
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P)'.  =  sin"1  { 


i  n 

•  -If  tu 


.  sin0  ) 
nr  (uu)  S' 


-  sin 


liq 

-1  f  1.4943 


IT6260  'sl"(48-8 


n;  =  43. 73° 


(4. 16; 


For  the  situation  that  k  is  parallel  to  pNLS(  one  has  e  =  48,  8° 


and  the  corresponding  n  can  be  computed  by  equation  (4.69 


2  !U 

n®  (0)  =  n_°  =  1.  5131 

2(ij  2(jj 


T 


as 


Thus  the  corresponding  angle  .f  incidence  fl"  is  given  by 


q"  =  sin*1  r  2,m 


I - 7—.  •  sinfl  [ 

(uj)  T> 

liq 


.  -1  f  1.  5131  „  o  ■) 

i  - Sln  tr6i^  • sm(48-8  )} 


=  44. 42° 
1 


(4.17; 


Due  to  the  intrinsic  birefringent  phenomena  in  KDP  crystal 

the  two  harmonic  beams  have  larger  spatial  s  eparation  when  q  become 

i 

large.  Furthermore,  information  from  equations  (4.  16)  and  (4.  17 )  indicate 
that  homonogeneous  and  inhomogeneous  harmonic  intensities  will  never 
be  simultaneously  zero  value  at  the  same  angle  of  incidence  9 

i 

Therefore  the  value  ITotal(2('))  given  by  equation  (4.14)  will  never 


be 


come 


,NLS 


zero  as  in  the  previous  case  where  the  nonlinear  polarization 


P  is  parallel  to  the  face  normal  direction,  i.  e.  ,  9  =0  =0.  The 

S  X 

experimental  points  are  indeed  in  good  agreement  to  the  theoretical 
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prediction  given  by  the  above  analysis.  It  is  noted  that  the  dynamical 
range  of  the  total  harmonic  intensity  for  this  case  is  about  10  orders  of 
magnit  ide  which  is  the  same  as  in  the  previous  case. 

The  investigation  relating  to  SHG  in  transmission  emanating  from  a  KDP 
crystal  immersed  in  an  optically  denser  linear  liquid  1-bromonaphthalene 
can  be  concluded  that  the  experimental  results  of  individual  as  well  as  the 
total  harmonic  intensity  are  in  good  agreement  with  the  theory  developed 
by  Bloembergen  and  Pershan  .  The  regime  of  transmitted  second 
harmonic  waves  corresponding  to  an  angular  range  from  normal  incidence 
to  critical  angle  '  as  been  demonstrated. 

c-  Second  Harmonic  Generation  (SHG)  by  Two  Beam  Spatial  Mixing 
(Noncollinear  Phase  Matched  Experiment). 

The  theory  developed  by  Bloembergen  and  Pershan(12)  has  also 

predicted  the  possibility  of  SHG  using  two  fundamental  beam  mixing 

inside  the  nonlinea  r  medium.  The  first  demonstration  of  SHG  by  two 

beam  spatial  mixing  (TRSM)  technique  was  performed  with  Q-switched 

ruby  laser  radiation  incident  on  GaAs  by  Ducuing  and  Bloembergen.  (25) 

(2  8 ) 

Later  Bloembergen,  Simon  and  Lee  performed  TBSM  experiment 

o 

using  Raman  laser  of  wavele  ngth  9770A  incident  upon  NaClC^  .  Those 
previous  experiments  were  performed  without  phase  matched  condition. 

In  the  present  study  of  TBSM  experiment,  an  uniaxial  KDP  crystal 
prepared  for  noncollinear  phase  matched  condition  was  employed  and 
the  fundamental  beam  was  Q-switched  pulse  from  the  Nd:glass  laser.  The 
purpose  is  to  extend  the  theory  into  the  domain  of  noncollinear  phase 
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matched  condition  from  an  uniaxial  KDP  crystal  and  as  a  consequence  the 
knowledge  of  polarization  and  SHG  by  TBSM  technique  will  be  utilized  for 

picosecond  pulsewidth  measurements  of  the  Nd:glass  laser. 

The  KDP  crystal  used  in  the  experiment  has  dimensions  of  25  x  15  x  8 

3 

mm  .  It  has  crystallographic  orientation  that  its  opfic  (z)  axis  is 
parallel  to  the  entrance  surface.  The  two  equal  intensity  fundamental 
beams  are  incident  from  opposite  sides  of  the  face  normal  which  is  in 
r  1 10]  direction.  The  angles  of  incidence  of  each  beam  are  equal  to  9 

i 

The  fundamental  beams  are  polarized  in  r  1  10]  direction  and 

-4 

according  to  equation  (2.31a)  the  nonlinea  r  polar ization  PNLSis  in  the  [  001  ] 
direction  (optic  axis).  The  polarization  and  crystallographic  orientation 
of  the  crystal  are  shown  in  the  inset  of  Fig.  4-10.  The  overall  experimental 
arrangement  and  the  transmitted  SHI  in  this  experiment  have  been  in 
Fig.  3-6  of  Chapter  HI. 

When  the  two  fundamental  beams  of  equal  intensity  impinge  on  the 

nonlinear  KDP  crystal  with  equal  ang’e  of  incidence  9.  from  the 

J  i 

opposite  sides  of  the  face  normal,  a  nonlinear  polarization  pNLS  is 
created  in  the  plane  of  incidence  with  zero  tangential  component  of  the  wave 
vector.  Since  the  polarization  has  the  same  phase  at  all  points  along  the 
surface,  it  radiates  second  harmonic  waves  in  the  normal  imvard  and 
outward  direction.  According  to  this  situation,  equation  (2.15)  can  be 
applied  directly  and  one  has 


eT  =  9S  -  9r  ■  0 


(4.18) 


The  value  of  g  inside  the  KDP  is  determined  by  the  length  of  the 


source  wave  vector  k  ; 

s 
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kc  ~  ■/  £_  a  -  n  cos9  a 

be  S  c  id  t 


(4.  19) 


t  . 


where  9  is  the  refracted  angle  for  the  individual  transmitted  fundamental 
b  ea  m . 


By  using  equation  (2. 30a)  and  n^  -  1.  0  ,  the  kg  given  by  equation 

(4.  19)  can  be  expressed  in  term  of  9  as 

i 


k  =  —  n°  ri-sin29./sin29Cr(uu)]1/2  a 
S  c  uj  -  ] 


(4.20) 


where  c r (fl) )  is  66.78°. 


According  to  equation  (4.  20),  kp  has  maximum  value  at  9=0°  and 

b  1 


c  r 

vanishes  when  9.  =  9  (uu)  .  This  is  in  agreement  to  the  fact  that  when 
the  two  fundamental  beams  are  incident  at  critical  angle  9Cr(uu)  the 
fundamental  beam  will  be  along  the  bourdary  of  the  interfaces  and  their 
normal  component  with  respect  to  the  interface  is  zero. 

•4 

From  equation  (4.20),  one  can  vary  the  length  of  k  which  is  along 

O 

-4 

the  face  normal  direction  such  that  it  is  equal  to  the  length  of  k^  in  the 
same  direction.  Thus  the  noncollinear  phase  matched  condition  is 
achieved. 

The  physical  description  which  described  this  condition  has  been  given 
in  chapter  II.  According  to  the  noncollinear  phase  matched  condition, 

one  has  for  normal  component  condition  given  by  equation  (2.51)  as 

ui  2ij u 

n  n 

'U  .  o  i  e  , 

n„  .  (jj  cosa  +  —  .  uucoso  =  -  (tt/2  )  .  2uu 
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n  coso  =  n  (tt/2 
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Thus 
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tu 


=  cos-1 


L  1.4943 


(4 


a'  =  10. 17° 


The  corres ponding  angle  of  incidence  or  the  noncollinear  phase 
matched  angle,  is  given  by  equation  (2.  53)  as 


m  _  .  -I-  0 


i  -  sm  n  sin^t' 
1 


-1  - 


-  sin  1.  4943  sin  10.  17°  1 


9,m  =  15.29° 

l 


The  total  transmitted  harmonic  intensity  IT<Jt (2tu)  was  observed 
in  the  direction  of  face  normal  of  the  exist  face  of  the  crystal.  It  was 
measured  as  a  function  of  angle  of  incidence  and  compared  to  the 
computed  theoretical  curve.  The  theoretical  curve  shown  in  Fig.  4-10 
was  computed  by  using  equation  (4.  15).  The  justification  of  using 


q  tion  (4.  15)  for  ITota](2t«l)  is  the  same  as  given  in  the  previous 


section.  The  anomalously  high  total  harmonic  intensity  I  (2w) 

Total 


m  Fig.  4-10  is  due  to  the  enhancement  of  nonlinear  Fresnel  factors 


NL  NL  m 

S’  11  and  FT,11  a‘L  °i  =  0i  •  This  can  be  mathematically  shown 


in  the  following. 

From  equation  (2.35b),  one  has 
NL  _  sin  ct 


S,  11  e  -  g 
eS  eT 


;  where  =  90c 


.  22  ) 
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(4.23) 


Furthermore,  since  9  =  9  =  9  =  0  and  a  =  90°  for  this  case 

o  1  K 

the  nonlinear  Fresnel  factor  F^^  given  by  equation  (2.35c)  is  reduced 


to 


,NL  _  _2S  1  _ 1_ 
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where  e  =  1.  0  for  air,  thus 
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Here  9  in  equations  (4.  22 )  and  (4.  23  )  is  the  same  pa  rameter  used 

in  equation  (4.19).  It  is  noted  that  both  F^^  and  F  NL  given  bv 

S,  11  T,  11  B 

equations  (4.  23)  and  (4.24)  re spectively  have  the  term  in  their  denominators 
satisfied  the  phase  matched  condition  given  by  equation  (4.12).  In  other 


NL 
S,  11 


,  NL 
T,  11 


words  F„  ,,  and  F.1’  ,,  become  infinite  at  9L  =  & '  =  10.  17"  .  Therefore 


upon  substituting  equation  (4.  23)  and  (4.  24)  into  equation  (4.  15)  one  has 
a  singularity  point  for  at  ^e  noncollinear  phase  matched 


condition. 


Note  the  striking  agreement  between  experimental  points  and  the 
computed  theoretical  curve  shown  in  Fig.  4-10.  The  noncollinear  phase 
matched  angle  9^  is  indeed  equal  to  15.29  as  predicted  by  the  theory. 
The  total  harmonic  intensity  IT  tal'2«j)  changes  by  about  four  orders  of 
magnitude  when  the  angle  of  incident  changes  by  a  few  tenth  of  a  degree 


from  9 


1 1 1 

.  .  This  is  for  the  first  time  that  the  theoretical  prediction  of 

the  enhancement  of  the  nonlinear  Fresnel  factors  F  and  F  ^ ^ 

S,  11  T,  11 

for  a  uniaxial  KDP  crystal  at  noncollinear  phase  matched  condition  has 
been  demonstrated. 


D.  Measurement  of  Picosecond  Pulsewidth 

The  agreement  of  experimental  data  with  the  theory  for  the  TBSM 
experiment  described  in  the  previous  section  encourages  us  to  utilize 
the  results  for  measurement  of  picosecond  pulsewidth  of  mode  locked 
Nd:glass  laser.  Since  KDP  is  a  nonlinear  crystal  from  which  phase 
matchable  SHG  can  be  produced,  then  it  is  very  convenient  to  use  its 
square  law  intensity  characteristic  to  measure  the  second  order  auto¬ 
correlation  function  of  Nd:glass  'aser  radiation. 

In  the  experiment  the  KDP  crystal  used  in  the  TBSM  experiment  in 
the  previous  section  served  as  the  nonlinear  medium.  The  KDP  crystal 
was  in  the  air  and  its  crystallographic  orientation  was  the  same  as 
indicated  in  the  inset  of  Fig.  4-10.  The  fundamental  beam  were  mode 


locked  pulse  trains  of  300-600  nanoseconds  duration  as  chown  in  Fig.  4-11 


and  their  polarizations  were  in  the  HlO  ]  di rection. The  two  fundamental  b 


earns 


were  incident  upon  the  crystal  at  noncollinear  phase  matched  angle  9  ,m. 
The  detail  of  experimental  setup  was  given  in  Fig.  3-6  of  Chapter  III. 
According  to  the  result  of  TBSM  experiment,  it  required  both  spatial  and 
temporal  overlapping  of  the  two  fundamental  beams  inside  the  crystal. 
The  temporal  overlapping  was  critical  and  could  be  achieved  only  by 
straight  forward  scanning  of  the  temporal  delay  of  one  beam  with  respect 
to  the  other  until  SHI  was  observed.  From  equation  (2.  62)  one  can 


see  that  the  transmitted  harmonic  intensity  I  (2tw )  is  directly 

1 0 13 1 

(2 ) 

proportional  to  the  second  order  auto-correlation  function  G  (t). 

(2 ) 

The  G  (t)  obtained  from  this  experiment  is  a  background  free  function 
since  a  single  fundamental  beam  cannot  generate  SHI  in  the  normal 
direction  to  the  exit  face  of  the  crystal.  In  the  experiment  the 

transmitted  harmonic  intensity  was  observed  as  a  function  of  temooral 

(2  ) 

delay  t  and  as  a  consequence  the  G  (t)  was  directly  mapped  out  as 
shown  in  Fig.  4-12  . 

The  resolution  in  the  picosecond  pulsewidth  measurement  is  stronglv 
dependent  on  the  thickness  of  the  KDP  crystal  used  in  the  experiment.  The  used 
crystal  thickness  of  8mm,  can  give  a  relatively  large  volume  of  interaction 
between  the  two  fundamental  beams  inside  the  crystal.  The  depth  of  the 
interacting  volume  inside  the  crystal  was  reduced  to  a  minimum  value  such 
that  the  SHI  could  be  comfortably  observed  for  the  entire  scanning  delay 
t  .  This  was  done  by  means  of  masking  the  entrance  surface  of  the 
crystal  by  a  0.  2  mm  wide  slit  ,  which  would  give  a  depth  of  interaction  of 
about  0.  6  mm  corresponding  to  a  time  duration  of  3.  0  picosecond  for  the 
light  travelling  inside  the  KDP  crystal.  From  the  experimental  result 
given  in  Fig.  4-12,  the  full  width  at  half  height  of  the  curve  G^\t)  gives 
about  2.67  picoseconds  for  the  width  of  the  auto-correlation  curve.  Since 
the  scanning  introduces  a  factor  of  two  for  the  temporal  delay  of  one  pulse 
in  relation  to  the  other,  thus  the  measured  pulsewidth  is  about  5  picoseconds. 

The  accuracy  of  this  measurement  can  be  improved  by  using  the  same 
crystallographic  orientation  of  a  KDP  crystal  with  a  thickness 
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of  the  order  of  only  0.1  -  0.2  mm.  Another  factor  that  contributes 
to  the  larger  pulsewidth  is  the  characteristic  of  the  laser  system  itself. 

It  is  an  intrinsic  problem  of  the  system  itself  whether  or  not  it  can 
reproduce  a  perfectly  mode  locked  pulse  train  .  Any  partially  mode 
locked  pulse  train,  which  cannot  be  distinguishable  by  the  fast  oscilloscope, 
will  yield  a  larger  measured  pulsewidth. 

Apart  from  the  resolution  limitation,  the  principle  of  noncollinear 
phase  matched  technique  via  square  law  of  SHG  can  be  used  for  second 
order  auto-correlation  mapping  and  for  picosecond  pulsewidth  measure¬ 
ment.  The  advantage  of  very  high  SHI  by  this  method  could  be  utilized 
in  the  future  for  the  measurement  of  the  pulsewidth  via  photographic 
methods  which  will  be  very  convenient  to  operate. 


CHAPTER  V 
DISCUSSION 

During  the  course  of  the  present  study  an  attempt  was  made  to  ex - 

(l2) 

tend  the  study  and  verification  of  the  theory  to  the  third  harmonic 
generation  (THG)  by  using  two  beam  spatial  mixing  (TBSM)  technique 
The  source  of  excitation  was  a  Q-switched  Nd:glass  laser  and  a  nonlinear 

optical  medium  immersed  in  an  isotropic  liquid  consisting  of  fuchsin  red 
dye  molecules  dissolved  in  hexafluoroacetone  sesquihydrate.  The  concen¬ 
tration  of  the  dye  solution  was  set  at  45  gm/liter.  At  this  concentration 
the  phase  matchable  THG  in  normal  incidence  and  total  reflection  of  a 
single  Ndrglass  laser  beam  was  demonstrated  earlier  by  Bey  et  al^3°’  5°\ 
and  indices  of  refraction  of  the  solution  were  given  as 

n(<jj)  =  n(3ou)  =  1.  3205  (5.1) 

The  experimental  setup  for  this  investigation  is  shown  in  Fig.  5-1. 

The  monitor  channel  was  similar  to  the  previous  setup  except  for  the 
fuchsin  pyrex  cell  and  an  interference  tilter  of  a  transmission  peak  at  3500  A° 
were  used  instead  of  the  quartz  platelet  and  second  harmonic  interference 
filter  respectively.  The  funda mental  beam  was  split  into  two  beams  of 
equal  intensity  by  the  beam  splitter.  Subsequently  the  two  beams  were 
incident  on  the  fuchsin  cell  with  equal  angle  of  incidence  from  the  opposite 
side  of  a  face  normal  of  the  cell.  The  signal  channel,  as  indicated  in 
Fig.  5-1  could  be  rotated  around  the  vertical  axis  which  was  tangential  to 
the  entrance  face  of  the  fuchsin  cell. 

The  signal  arm  was  scanned  through  a  large  interval  of  angle  by  1° 
for  each  step.  Despite  careful  alignment  of  every  optical  component, 
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varying  the  angle  of  incidence  9.  and  repeated  scannings,  the  transmitted 

third  harmonic  intensity  (THI)  was  not  observed.  The  reasons  that  justify 

the  result  can  be  mathematically  explained  in  terms  of  a  noncoliinear  phase 

matched  picture  which  will  be  given  belcw. 

The  two  incident  beams  of  equal  intensity  are  incident  from  quartz 
12 

with  3.-0.  -  9.  as  shown  in  Fig.  5-2.  Since  the  dye  fuchsin  solution 

is  an  isotropic  medium, the  locus  of  transmitted  third  harmonic  wave 
vector  kT  will  be  a  semicircle  of  radius  3  |  kfc  |  inside  the  medium. 

-4 

The  kfc  is  the  transmitted  fundamental  wave  vector.  According  to  the 
phase  matched  condition  given  by  several  author s  3 '  34,  43  \  the  three 
transmitted  fundamental  wave  vectors  must  add  up  so  that  the  resultant 
vector  is  equal  to  k^w).  Here  one  assumes  that  one  fundamental  beam 
contributes  the  transmitted  fundamental  wave  vector  2k  (yj)  whereas 

-4 

the  other  is  k^(uj)  .  Furthermore,  boundary  conditions  given  by 
Bloembergen  and  Pershan  require  that  the  tangential  components  of 
those  wave  vectores  must  be  conserved.  By  using  equation  (2.  49)  one 
obtains  from  Fig.  2-5  that 
tangential  component: 

T 2 k^_ (uj )  -  k^fuj)]  sin9^  =  k^,(3uj)  sin9,j,  (5.2) 

normal  component: 

[2k^(uj)  +  k^(uj)]  cos9^  =  k^(3uu  )  cos8^,  (5_  3) 

Since  k_(u) )  =  -  n(w)  and  k^uj)  =  —  n(3uu)  ,  equation  (5.  2 )  and 

(5.  3)  can  be  deduced  to  equation  (5.  4)  and  (5.  5)  respectively  as 
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(5.4) 
(5.  5) 


Since  the  dye  fuchsin  solution  has  a  concentration  such  that 

n(jj)  =  n ( 3 ix) )  ,  the  solution  of  the  simultaneous  equations  (5.4)  and  (5.  5) 
becomes 


(5.6) 


This  is  not  a  surprising  result.  It  has  the  physical  meaning  that 

phase  matched  THG  will  be  created  if  the  fundamental  beam  is  normally 

im  ident  upon  the  dye  fuchsin  solution.  Under  this  situation  one  obtains 

rom  equation  (5.  3)  that  3k^  =  and  the  two  fundamental  beams  will 

collapse  into  a  single  beam.  This  situation  has  already  been  investigated 
(50) 

by  Bey  et  al 

From  the  above  analysis  it  is  noted  that  when  the  dye  fuchsin  of 
concentration  45  gm/liter  is  employed  as  a  nonlinear  medium,  the 
phase  matched  THG  is  produced  under  normal  incidence  of  the  fundamental 
beam.  However,  in  the  experimental  investigation  the  angle  of  incidence 

was  varied  in  2°  steps  from  10  to  25  .  Thi s  variation  of  the  angle  of 

incidence  apart  from  0  will  yield  phase  mismatched  condition  and  since 

the  dye  fuchsin  solution  is  strongly  absorbed  at  the  third  harmonic  line 
o 

(3530A),  a  transmitted  third  harmonic  intensity  will  not  be  observed 
in  the  experiment. 

From  equation  (5.4)  and  (5.  5),  future  analysis  can  be  made  if  one 
selects  a  dye  concentration  different  from  45  gm/liter.  One  has  to 
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use  higher  dye  concentration  due  to  the  fact  that  n(uu)  is  greater  than 
n(3'.jj )  in  this  region  and  the  equation  (5.4)  and  (5.  5)  will  give  physical 
solution.  This  can  be  mathematically  demonstrated  in  the  following 
example. 

Ry  using  the  dye  concentration  of  50gm/liter  the  indices  of 
refraction  at  fundamental  and  third  harmonic  frequencies  will  be  given 
as 


n(-t')  =  I.  3225 


n(3;ij)  =1.3210 


From  equation  (5.  4)  and  (5.  5)  one  obtains 

i  =  (9Si„T  +  cos2, 

T  T  1  n(uu) 


,  ,o  ■  2n  ,  ,,  f  1.  3210 
(8sin  9t+  1)  1T732Z5 


Q  =  0.  97 
T 


(5.7) 


(5.  8) 


Thus,  from  equation  (5.4)  the  corresponding  angle  of  incidence  in 
Fig.  5-2  will  be 

P.m  =  2.88°  (5.9) 

This  is  the  phase  matched  angle  for  a  dye  concentration  of  50  gm/liter. 
One  can  use  even  higher  dye  concentration  in  order  to  achieve  a  relatively 
large  phase  matched  angle.  However,  due  to  a  lack  of  knowledge  of  the 
indices  of  refraction  at  hand,  the  experimental  investigation  has  not  been 
carried  out.  The  above  analysis  will  serve  for  future  experiment  involving 
TBSM  for  third  harmonic  generation  from  isotropic  nonlinear  media 


with  n(ou  )  >  n  ( 3(ju  )  . 
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Even  though  the  experimental  result  for  the  case  of  dye  concentration 

of  45  gm/liter,  turned  out  to  be  negative, it  can  be  explained  by  the  theo 
,  ,  ,  ,  .  .  (12,  33,  34,  43 ) 

developed  by  several  authors  In  a  relatively -new  field  it 
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cannot  be  expected  that  experimental  results  always  agree  with 
theoretical  predictions  ,  However,  to  be  able  to  ultilize  theoretical 
knowledges  to  explain  the  phenomena  observed  in  the  experimental 
investigation  will  result  a  better  understanding  and  direct  contribution 
to  the  theory. 

For  this  reason  it  is  very  important  that  experimental  and  theoretical 
go  hand  in  hand  for  increase  our  knowledge  of  <:he  field 


CHAPTER  VI 
CONCLUSION 


The  nonlinear  optical  interaction  in  a  KDP  prism  in  relation  to 
second  harmonic  generation  was  investigated  by  means  of  Q-switched 
and  mode-locked  Nd:glass  laser  pulses  and  the  results  obtained  from 
the  investigation  were  subsequently  utilized  to  measure  the  second  order 
intensity  auto-correlation  function.  The  experimental  results  obtained  from 
this  investigation  provide  additional  rather  striking  and  detailed  confirma¬ 
tion  of  the  theory  developed  by  Bloembergen  and  Pershan(12)  on  the 
behavior  of  the  light  waves  at  the  boundary  of  a  nonlinear  crystalline  KDP 
medium.  Although  the  theory  was  developed  for  a  homogeneous  incident 
plane  wave  of  infinite  cross  section,  the  results  are  well  described  by 
taking  the  field  solution  for  the  infinite  plane  wave  and  the  cutting  off 
the  beam  by  modification  of  the  aperture  width  expression  and  taking 
linear  Fresnel  factor  into  consideration.  The  recasted  theoretical 
expressions  given  in  chapter  II  have  been  well  verified  and  are  valid  for  an 
uniaxial  crystal,  e.  g.  ,  KDP.  The  reasons  that  the  modification  of  the 
infinite  plane  wave  case  works  so  well  is  due  to  the  fact  that  the  dimension 
of  the  relevant  experimental  parameters  are  large  compared  to  the 

characteristic  dimension  of  the  physical  problem.  Detailed  justification 
has  been  given  by  Simon^^. 

The  results  obtained  from  the  experiment  of  phase  matched  second 

harmonic  generation  at  total  reflection  are  in  agreement  with  the  other 

,  ..  (88) 

observations  under  the  same  condition  except  different  wave  length 
and  type  of  the  fundamental  beam.  When  the  KDP  crystal  being  phase 
matched  along  its  surface  was  rotated  about  its  face  normal  by  180°  it 
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was  found  that  the  reflected  second  harmonic  intensity  in  the  neighborhood 
of  critical  angle  was  always  tower  .  This  phenomena  can  be  qualitatively 
explained  by  an  intrinsic  property,  the  so  called  double  refraction  which  re¬ 
sults  m  a  walk-off  effect  in  the  uniaxial  crystal.  However,  an  exact 
mathematical  treatment  that  accounts  for  a  detailed  explanation  of 
this  phenomena  remains  for  future  investigations. 

With  good  utilization  of  the  advantages  of  mode  locked  Nd:glass  laser 
pulses  for  the  observation  of  reflected  second  harmonic  intensity  fr6m 
KDP  crystal,  the  Nonlinear  Brewster'  s  angle  of  the  transparent 
medium  (KDP)  has  been  first  demonstrated  with  an  excellent  agreement 
to  the  theoretical  prediction. 

For  the  transmitted  second  harmonic  waves,  the  angular  dependence 
of  homogeneous,  inhomogeneous  as  well  as  the  total  harmonic  intensities 
emanated  from  KDP  crystals  were  experimentally  confirmed  with  the 
theory.  The  phenomena  of  phase  matched  SHG  along  the  oblique  direction 
by  means  of  birefringence'33'  34)  was  investigated  and  the  results  were  also 
in  good  agreement  to  the  theory.  The  relative  differences  of  SHI  between 
the  maximum  and  minimum  as  indicated  in  Fig.  4-8  and  4-9,  were 
observed  to  be  about  ten  orders  of  magnitude.  The  striking  agreement 
of  the  experimental  result  of  noncollinear  phase  matched  SHG  by  means 
of  two  beam  spatial  mixing  (TBSM)  to  the  theoretical  prediction  was 
utilized  to  measure  the  second  order  intensity  auto-correlation  function 
of  Ndrglass  laser  radiation.  This?in  principle,  can  be  used  for  the  construc¬ 
tion  of  a  device  for  picosecond  pulsewidth  measurement  using  simple 
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photographic  techniques.  This  is  due  to  the  fact  that  the  SHI  emanating 

from  the  KDP  crystal  under  noncollinear  phase  matched  condition  is 
so  intense  that  it  can  be  photographed. 

Finally,  the  results  obtained  from  this  investigation  confirm  and 
extend  the  theory  of  Bloembergen  and  Pershan(12)  for  -onlinear  optical 
interaction  in  an  uniaxial  crystal  (KDP)  in  the  regimes  of  reflected  and 
transmitted  second  harmonic  intensities. 


APPENDIX 

COMPUTER  PROGRAMS  I  OR  THEORETICAL  CURVES 
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■KUN  /B.3U5168.BHANTHUMNA, 1.100 
•PASSWD  17B9FI 
•  FOR  ♦  I S  MAIN 

COMPLEX  Alt  A2,  A3.  A4,  A5.  A6.  A7.  A8 .  FNL »  FLT  *  C.  D 

C  r  (  1 . 0 , 0 . 0  ) 

0  =  (J. 0.1.0) 

READ  ( 5 • 7 ) A 
7  FORMAT ( F  8 . 5  ) 

10  XINC  =  0.1 
A  =  A+XINC 
PI  =  3.1416 

IF  (  A-75.0  I  12,  12,  100 

12  AX  =  A*P I / 1 80  •  0 
Z  =  SIN(AX) 

A 1  =  Z«C 
Z  =  COS ( AX ) 

A2  =  Z*C 
B1  =  1.513122 
B2  =  1.47122 
B3  =  1.494305 
01  =  1.6260 
02  =  1.670 
Z  =  (01/02 ) * A 1 
A3  =  Z*C 

Z  =  SORT  (  1 . 0—  (A3)  ■** 2 •  0 ) 

A4  =  Z*C 

13  IF  (A-64,797)  14,  14,  15 

14  Z  =  (  D1/B3)  *A 1 
A5  =  Z*C 

Z  =  SORT  (  1.0-  (A5)**2.0) 

A6  =  Z#C 
Z  =  ()1/B2)*A1 
A 7  =  Z*C 

Z  =  SORT  (  1.0-( A7)**2.0) 

A8  =  Z*C 
GO  TO  80 

15  IF  (A-66.781)  16,  16,  17 

16  Z  =  <D1/R3)*A1 
A5  =  Z*C 

Z  »  SORT  (  1.0-  <A5)**2.0) 

A6  =  Z»C 
Z  =  (D1/B2I *A1 
A 7  =  Z#C 

Z  =  SORT  ( (A7)**2, 0-1.0) 

A8  =  Z*D 


GO  TO  80 

17  Z  =  (01/B3 I *A1 
A5  =  Z*C 

Z  •=  SORTt (A5>**2. 0-1.0) 

A6  *  Z*D 
Z  =  (D1/B2)*A1 
A 7  =  Z*C 

Z  =  SQRT( (A7)**2. 0-1,0) 

A8  =  Z*D 

80  FLT  =  ( ( 2.0) *A  2*  A  5  )  /  ( A  1*  A6  +A2*A5) 

FNL  =  (A5«A7*A7«A7)/(A3*(A7*A4  +A8*A3)*(A8*A4  +  A7*A3) 

1*<A7*A6  +  A8*A5  )  ) 

FLTR  =  REAL ( FLT ) 

FLT I  =  AIMAG(FLT) 

FNLR  *  REAL  (FNL) 

FNLI  =»  AIMAG  (FNL) 

FNLM  =  SORT  (  F f(LR *FNLR  +  FNLI*FNLI) 

FLTM  ^  SORTt  FLTR*FLTR  +  FlTI#FLTI) 

SHI  =  (FNLM**2.0)*(FlTM**4.0)#(A4/A2) 

WRITE  (6,90)  A,  SHI 

90  FORMAT  (1H0,  •  INC.  ANGLE  =  *,F6.3»  '  REFLECTED  SHI  =  *F10 

GO  TO  10 
100  END 
•  XOT 
59,9000 
•FIN 


Computer  Program  for  the  Theoretical  Curve 
of  Figure  4-1 
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i 

I 

I 

l 
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I 
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•RUN  V8 • 305 16802 • BHANTHUMNA * 3*200 
•FOR. IS  MAIN 

COMPLEX  F» A1 .A2.A3 .A4.A5.A6.F1 .FNL.FL.C.D 
C*I1. 0,0.0) 

0=10,0,1,0) 

IFIA-/5.0)  12.12,50 

12  IFIA-66.615)  13,13.14 

13  X=SIN( A*P1/180.0) 

X  =  1 . 089  2*X 

X  =  X*X 

X*S0RT< 1,0-X) 

F  =  X*C 
GO  TC  15 

14  X=SIN< A*PI /180.0) 

X=  1.0892*X 
X  =  X*X 

X*SORT ( X-l .0 ) 

F*X*0 

15  AX  SA*P1 / 1 80,0 
2»SINIAX) 

Fl=  1 «  0892#Z*C 
A1»2*C 
ZnCOS ( AX ) 

A2*Z»C 

ALPHA“41,2*P 1/180.0 
Z»SIN< ALPHA) 

A3  =  Z*C 
Z«COS( ALPHA) 

A4«=Z*C 

R*  AS  IN ( 0. 9746*SI N ( Ax ) ) 

Z  =  S IN ( R  I 

A5=Z*C 

2»COS(R) 

A6=Z*C 

FNL=(F1*F1*F1I*(2.0*F1*F*A4-A3+2.0*F1*F1*A3)/IA5*I2.0*F1*F)*(F1*A6 
1+F*A5)«(F*A6+F1*A5) ) 

FL=2.0*Fl*A2/( A1*F+A2*F1) 

FNLRsREALIFNL) 

FNLI=AIMAG(FNL) 

FLR=REAL<FL) 

FLI=AIMAG<FL) 

FNLM=50Rt«FNLR«FNLR+FNLI*FNLI ) 

FLM  =  SQRT(FLR*FLR+FLI*FLI  ) 

SHI=(FNLM**2.0)*(FLM**4.0) 

WRITE16.20)A.SHI 

20  FORMAT ( 1H0. *  INCIDENT  ANGLE  «»F8.5.'  INTENSITY  OF  SHG  =»F10.3) 

50  FND 
59.90000 
00 


J 

I 


Computer  Program  for  the  Theoretical  Curve 
of  Figure  4-2 
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•RUN  VB.30516802»BHANThUMNa*3*200 
•FOR* IS  MAIN 

COMPLEX  A1.A2.A3,A4,a5»A6.A7,a8,FNL.FL.C,D 
C*  (1. 0*0.0) 

D=  10.0*1.0) 

READ ( 5  *  7 ) A 
7  FORMAT ( F8.5 ) 

10  XINC  *C.l 
A  =  A+XINC 
PI  =  3. 1416 
IF  (A-75.0)  12.12*100 

12  AX=  A*P I / 1 80 • 0 
Z=  SIN(AX) 

Al=  Z*C 

Z  =  COS(AX) 

A2=  2*C 

Z=  (1. 6277/1. 670)*A1 
A3=  Z*C 

2=  SORT! 1.0- I (1.62 77/ 1.670 )*A1 ) **2.0) 

A4  =  Z*C 

IF  IA-66.61)  13.13.U 

13  Z=(l. 5277/1. 4942)*A1 
A5=  Z*C 

2=  SC  7T ( 1 • 0- ( ( (  1, 6277/1. 49a2)*A1)** 2.0) ) 
A6=  Z*C 
GO  TO  15 

14  Z=  ( 1.6277/1. 4942)*A1 
A5  =  Z*C 


Z*=  SORT  (  (  (  (  1.6277/ 1.4942) *A1  >**2. 01-1.0) 

A6=  Z*D 

15  ALPHA3  48. 8*P 1/180.0 
Z=  SI N( ALPHA ) 

A7=  Z*C 
Z=  COS( ALPHA) 

A8a  Z*C 

FNL=  (  (  A5*A5*A5)*(  A5*A8+A6*A7)  )  / ( A3* I  A5*A6+A  *,*A5  ) * (  A5* A4+A6*A3  ) 
1*(A6*A4+A5*A3)) 

FL=  2.0*A5*A2/(A1*A6+A2*A5) 

FNLR=  REAL(FNL) 

FNL I =  AIMAG(FNL) 

FLR=  REAL(FL) 

FL I =  AIMAG(FL) 

FNLM=  SORT(FNLR*FNLR  +  FNlI*FNLI) 

FLM=  SQRT(FLR*FLR  +  FLI*FU ) 

SHI  =  (FNLM**2.0)*(FLM**4,0 ) 

WRITE  (6*20)  A,  SHI 

20  FORMATIlHO.'  INCIDENT  ANGLE  *>*F8.5»'  INTENSITY  OF  SHG  «»F10.3) 
GO  TO  10 
100  END 
'  XOT 

59.90000 

•FIN 

•FIN 


Computer  Program  for  the  Theoretical  Curve 
of  Figure  4-3 


RUN  VS>3C516802.nHANTHUMNA.3.200 

'  FOR  >18  f’A  I  N 

READ ( 3  >  7 ) A 
7  FORMAT ( F8.S  ) 

10  XINC  =0.1 
A  =  A+X1NC 
PI  =  3.1416 
IF  ( A  —  50.0  I  12.12.100 
12  A X=  A*PI/180.0 
H  =  SIN ( AX  1 
A 1  =  2 
1=  COS ( AX ) 

A  2  =  2 

Z= ( 1.6277/1.6701*41 
A3=  2 

Z  =  SQRT (  1.0-1 ( 1. 6277/1. 6701*A1 1**2. 0) 

A4=  2 

Z-  (  1 .62 77/ 1 .4942 1 *A1 
AS  =  2 

2  -  SORT  (1.0-(( (1.6277/1.49421*411**2.0) 1 

A6=  2 

ALPHA  =  48.P*PI/180.0 
Z-  SIN(ALPKA) 

A  7  =  2 

Z=  COS(ALPHA) 

A8=  2 

2=  (1.6277/1.47161*41 
A9  =  2 

Z  =  SQR  T ( 1 ,0-( (  ( 1.6277/1  .47161*41  1**2. 0  1  1 
A  10=  2 

l*?Mi*A4^J*X3ir*Aft’Al0#A7”/(A3*(A9*A6  +  A10#A5l#'A9*ADDA10*A3> 

FL>  <2.G*A  2*  \S  1  /  (  A  1  *  «t>  +  A2*  AS  I 
SHI  =  ( FNL*«  1.0 1* ( Fl»*4.0 1 
WRITE  (6.20)  A.  SHI 

20  GO*™Tir°’’  I,NCI0ENT  ANGLt  =  'F'8-5’’  INTENSITY  OF  SHG  •  '  F 1 0 . 3  1 

100  END 
'  XOT 

34.90000 

'FIN 


Computer  Program  for  the  Theoretical  Curve 
of  Figure  4-6 


4> 


V 


•RUN  VB. 305168, BHANTHIJMNA, 1,100 
•PASSWD  17B9FI 
'FOR  *  IS  MAD 

READ ( 5 , 7 ) A 
7  FORMA  f ( F8,5  ) 

10  XINC  a  0.1 
A  *  A+XINC 
PI  =  3.1416 
IF  (A-90.0  )  12,12,100 
12  AX  =  .\*P  J  / 180.0 
A 1  =  SIN(AX) 

A2  =  COS ( AX  ) 

B2  =  1 ,  A  7 1 2  2 
B 1  =  1.513122 
B3  =  1.494305 
01  =  1.6260 
02  *  1.670 
A3  =  (B1*A1J**2.0 
A4  =  ( B2*A  2 ) **2.0 
A5  *  SORT ( A4+A3 ) 

B4  =  (Bl*B2)/A5 
CX  =  AX 
Cl  =  SIN(CX) 

C2  =  ASINI(B4*C1I/B3 J 


C3  *  <C2*180.0)/PI 
C4  *=  ASINI(B4*C1)/01) 

C5  =  <C4*180.0! /PI 
Cfc  =  ASIN( (B4*Cl)/02) 

C7  «  ( C6*l 80 , 0 ) /P I 
R1  =  SIN ( C4 ) 

R2  =  C0S(C4) 

R3  =  SINIC2I 
R4  *  COS(C2 I 
R5  =  3INIC6) 

FLT  =  (2.0*R2*R3)/(Ri*R4+R2*R3) 

FNLS  a  ( (R3)*C-1.0) >/( (B3) **2.0- (64) **2.0) 

™LT  *  (  (1.0)/ID2*B4j  )*IR3*A1*A1*A1)/(R5*(SI  I(  AX+C6  I  I  *  (  COS  (  AX-C6  )  * 
1#(5IN( AX+C2 1 1  I  -  !FNlS)*(B3/B4! 

SHI1  ,  B3#(FLT**4.0)*(FNLS**2.0)*R4/R2 
SHI  2  *  B4*(FLT**4.0)*(FNlT**2.0)*A2/R2 

SHI  =  SHI  1  +  SHI  2 

WR I TE I  6 , 90 )  A,C3,C5,  SHI  1 ,  SHI  2  ,  SHI,B4 


ZETAT  =  •  ,F6. 3  ,  •  ZETAS 


» F6, 3 i 


90  FORMAT  (1H0,  • 

1'  INC«ANGLE=  *  ,F6  •  3 , '  SHI1-  *,F10.3,*  SHI2«  *^10^3, 
1*  SHI  =  • ,F10.3 ,  »  INDEX  E«  '.F10.7I 
GO  TO  10 
10{  END 

•  XOT 
0»00f 00 

•  FIN 


Computer  Program  for  the  Theoretical  Curve 
of  Figure  4-7  and  4-8 


•RUN  VB. 305168, BHANTHUMNA, 1,100 

*  PASSWD  17B7FI 
•FOR,  IS  MAPI 

READ { 5 , 7 ) A 
7  FORMA  T ( F8.5  ) 

10  XINC  «  0,1 
A  «=  A+XINC 
PI  *  3.1416 
IF  (A-66, 81 12,12,100 
12  AX  =  .\*PI/180.0 
A1  =  SIN(AX) 

A2  •  COS( AX) 

B3  *  1.494305 
B4  *  1.47125 
01  =  1.6260 
A3  «  ASINI <D1*A1)/B3) 

A4  =  SI N  <  A3 ) 

A5  =  f:0S(A3) 

FLT  =  (2.0*A2*A4)/(A1*A5  +A2*A4> 

FNLS  *  ( ( A4  I  *  (-1.01 1 / ( (B3) **2.0  -(B4)**2.0) 

SHI1  »  B3*(FLT**4.0)*(FNtS**2.0)*A5/A2 
WRITE  (6,90)  A,  A3 ,  SHI  1 

90  FORMAT  (1HO,  •  INC. ANGLE  *  *,F6.3,  •  2ETAS  «  >,F6.3, 

1‘  SHI1  =  SF10.3) 

GO  TO  10 
100  END 

•  XQT 

64.50000 


Computer  Program  for  the  Theoretical  Curve 
of  Figure  4  -7  (  partial  ) 
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I 

I 

I 

f 

I 

I 

1 


i 


I 

I 


•RUN  VB.30516B.BHANTHUMNA, 1.100 
•PASSWD  1 7RJF 1 
•FOR, IS  MAIN 

READ! j, 7) A 
7  FORMAT ( F 1 2 • 8 ) 

01  •  1,4889679 
02  ■  1.4929732 
03  *  1.4968934 
04  «  1.5041322 
05  *  1.5097313 
06  •  1.5127767 
07  *  1.5131219 
08  •  1.5125449 
10  XINC  •  0.1 
A  «  A+X I NC 
PI  •  3.1416 
B1  *  1.6260 
B2  «  1.494305 
B3  *  1.670 
B4  «  (48. 8*P1 1/180.0 
IF  (A-5.0112.12.15 
12  YINC  *  0.000079328 
01  ■  01+yINC 
OX  ■  01 
60  TO  80 

15  IF  (A-10.01  16,16.25 

16  YINC  -  0.00007878 
02  ■  02  +  YINC 

OX  -  02 
GO  TO  80 

25  IF  (A-20.0)  26,26.35 

26  YINC  «  0.00007256 
03  «  J3+YINC 

OX  -  03 
GO  TO  80 

35  IF  (A-30.0)  36.36.45 

36  YINC  ’  0,00005881 
04  .  J4+YINC 

OX  *  04 
GO  TO  80 

45  IF  (A-  40.0)  46.46,55 

46  YINC  *  0.00003050 
05  *  J5+YINC 

OX  •  05 
GO  TO  80 

55  IF  (A-44.3)  56.56.  65 

56  YINC  •  0.00000793 
06  '  D6+YINC 

OX  -  06 
GO  TO  80 

65  IF  (A-50.0)  66.66.75 

66  YINC  -  0.00001034 
07  ■  07-YINC 

OX  ■  07 
GO  TO  80 

75  IF  (A-60.0)  76.76,100 

76  YINC  *  0.00004607 
08  ■  D8-YINC 

OX  «  08 

80  AX  *  A*°I / 1 80.0 
A 1  *  SIN(AX) 

A2  *  COS ( AX ) 

A3  «  <B1*A1)/B2 

A4  *  SORT  < 1,0-1 A3*A3) ) 

A5  *  (  B 1 • A 1 )  /DX 

A6  •  SORT ( 1.0-1 A5«A5) ) 

A 7  ■  (B1*A1)/B3 

A8  '  ASINI (B1*A1 ) /B2) 

A9  ■  AS1 N  < (B1*A1 ) /DX) 

A10  «  ASlN< <B1*A1 )/B3> 

FLT  ■  (2.0*A2*A3)/(Al*A4  ♦  A2«A3) 

FNLS  ■  (  SINIB4-A8)  )/(  (  B2 )  **2 . 0- (DX  . 0 ) 

FNLT  «  (  1.0/(81‘*DX))*(A3*A5*A5*SIN(A9*b4)  )  /  <  A7*S!N(  A9+A10  ) 

1*C0S(  A9— A10 )*SIN( A8  +  A9) )  - ( B2 /DX ) * ( S IN ( B4-A8 ) ) / ( (B2 )**2.0- 

1 ( DX ) **2,0  ) 

SHI  1  -  B1*(FLT**4.0)»(FNLS»*2.0)«A4/A2 
SHI  2  •  0X*<FLT«»4.0)«(FNLT**2.0)*A6/A2 
SHI  «  SHI1  +  SHI  2 

WRITE  (6,90)  A  «DX ,  SHll.  SHl2,  SHI 

90  FORMAT  <1H0,  '  INC.  ANGLE  ■  '.F6.3.  •  I.1DEX  ERAY  ■  »,F10.7. 

1*  SHU  ■  '  .FI  2,3,  •  SHl  /  ■  •  ,F12 . 3,  •  SHI  «  '.F12.31 

GO  TO  10 
IOC  END 

•  XQT 

0.00C300 

•  FIN 
•FIN 


4 


Computer  Program  for  the  Theoretical  Curve 
of  Figure  4-9 


•RUN  V8.30516802* BHANTHUMNA *  3  *  200 
•FOR* IS  MAIN 

READ ( 5  * 7 ) A 
7  FORMAT  (F8.5) 

10  XINC  =  0.1 
A  =  A+XINC 
PI  =  3,1416 
IFIA-25.0) 12.1 '.100 
12  AX=A*P 1/180.0 
Z=SINI AX) 

A1  =  Z 

Z  =  COS  t AX ) 

A2  =  Z 

Z=( 1.0/1. 4942)*A1 
A3  =  Z 

Z  =  SOR  T ( 1 , 0-A3*A3 ) 

A4  =  Z 

FNLT  =< 1.494 2/1.4712 )*A4*( 1.0/ ( ( 1 .4942* A4 ) **2 .0- ( 1 . 471 2 ) **2 . 0 ) ) 

1 -Ml. 0/1. 4712) 

FL  T*=  I2.0*A2*A3»/(A1*A4  +  A2*A3) 

FNLS  =  ( 1 . 0/ ( ( 1.4942 *A4)**^,0-( 1,47121**2.0) ) 

SHI1  =  (FLT**4.0)*(FNLT**2.0)/A2 
SH I  2  *  (FLT**4.0)*(FNLS**2.0)/A2 
SSH I  =  SH 1 1  +  SH 1 2 
WR I TE ( 6 • 20 )  A,  SSH1,  SHI  1  .  SHI2 

20  FORMAT ( 1H0, ’  INCIDENT  ANGLE  ='F8.5.  'TOTAL  INTENSITY  OF  SHG  =  ' 
IF  10. 3  »  • SH 1 1  =  'Flo. 3,  ' SH I  2  =  '  F10.3) 

GO  TO  10 
100  END 
•  XQT 
0.00000 
•FIN 


Computer  Program  for  the  Theoretical  Curve 
of  Figure  4-10 
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